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ABSTRACT 


In direct-sequence spread-spectrum systems, successful communications require 
phase synchronization of the m-sequence in the incoming signal with a locally generated 
m-sequence at the receiver. Many acquisition schemes which extract the phase of an 
incoming m-sequence have been studied, but most of them assume coherent demodu- 
lation (which is usually not available during acquisition) and/or independent samples 
(which introduce a loss in the effective signal to noise ratio (SNR)). 

This thesis investigates the performance of two m-sequence acquisition schemes in 
the presence of fading and data modulation. A fixed sample size test and a truncated 
sequential test are studied without the usual assumptions of coherent demodulation or 
independent samples. The effects of fading and data modulation on our schemes’ 
probability of false alarm, probability of detection and test length are thoroughly ex- 
plored. We find that channel fading in effect induces a loss of signal SNR, but the de- 
sired power of the tests can be restored by suitable adjustments in the decision processor. 
We find that the effects of data modulation are less severe, but more problematic to 


correct. 
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1. INTRODUCTION 


This thesis explores the effects of fading and data modulation on sequential 
pseudonoise (PN) acquisition schemes. We begin by providing a brief introductory 
overview of spread spectrum communication systems. The purpose of this introductory 
chapter is twofold. The inclusion of this material enhances the degree to which this 
thesis is self-contained. Additionally, the notation introduced in this chapter will be used 
in subsequent sections. 

A complete presentation of the theory and applications of spread spectrum com- 
munication systems would fill several volumes. The intent here is neither to be com- 
prehensive in scope nor specific in detail. We present, in this chapter, only a thumbnail 
sketch of the subject, including only those details necessary to establish a foundation for 
the following chapters. 


Spread spectrum signals allow us to [Ref. 1 : pp. 800-801]: 


l. overcome or minimize the effects of interference, whether it be intentional (e.g., 
jamming), unintentional (e.g., interference from other channel users) or self-induced 
(e.g., multipath). 


2. hide a signal by “burving” it below the background noise level, thus concealing the 
transmission from unintended listeners and eavesdroppers. 


3. achieve message privacy. 


In a conventional communication system employing amplitude modulation (AM) 
or phase modulation (PM), the communication signal bandwidth is approximately equal 
to twice the rate of data transmission. In spread spectrum (SS) systems the bandwidth 
of the communication signal is much greater than the data rate. In all SS systems, ran- 
domness is used to make interception and copying more formidable. 

Prior to 1980, the field of spread spectrum communications fell almost entirely 
within the military domain. Techniques were developed to allow communications which 
were not detectable by adversaries and to combat the disruptive effects of jamming by 
an opponent. In 1983 the Federal Communication Commission (FCC) opened, for the 
first time, three frequency bands for commercial spread spectrum use [Ref. 2]. Cellular 
and indoor commercial networks using spread spectrum techniques are now becoming 
a reality. 


In the following we will discuss several typical spread spectrum techniques. 


A. TIME HOPPING (TH) 

In a time hopping scheme, we transmit the communication signal at randomly cho- 
sen times for random burst durations. ЈЕ 7 is the fraction of time the transmitter is active 
(у < 1) and r, is the steady data rate from the user, then r, , the transmission rate when 


the transmitter 1s activated, is given by 
Po 
Р = | 


Note that the bandwidth of the signal that gets transmitted is proportional to r, (which 
may be much larger than r, ) so we have distributed, or spread, the signal over a wider 
frequency band. 

The primary advantage of TH is an improvement in the Signal to Jamming Ratio 
(SJR). Let S and J be the average signal power and average jamming power, respec- 


tively, at the receiver input. Then the SJR for a non-SS system is given by 
S 
5ЈЕ == (1,2) 


For a TH svstem, let S be the input power at the receiver when a burst transmission 


occurs. Then $= Sn and the SJR during burst reception 15 
зук = (Ss am 


Thus we see that the SJR in a TH svstem is better than that for a non-SS system. The 
amplification factor (= =) is called the processing gain (PG). 
The improvement in SJR in a time hopping system 1s countered by several disad- 
vantages: 
1. a data buffer must be employed in both the transmitter and receiver 
2. the TH transmitter must radiate at higher power levels 


3. the TH receiver must resynchronize on each new burst 


B. FREQUENCY HOPPING (FH) 

As previously noted, in all SS systems randomness is used to make the communi- 
cations system more immune to signal detection and jamming. In a system employing 
frequency hopping, the transmitter carrier frequency changes value (hops) as a function 


of time. If we let r, be the amount of time the carrier remains at each of the hopping 


frequencies, then the bandwidth of the carrier is always constant and equal to 2/1, = b 
Hz. If Bis the total range of frequencies we are assigned to hop over, then the number 
of distinct carriers used 1s equal to B/b. 

A FH receiver despreads the received signal by multiplying it by a locally generated 
hopping signal. Note that the hop sequence of the receiver’s local oscillator must match 
the hop sequence of the transmitter in order to recover the signal. 

As with time hopping, a major advantage of FH systems is an improved immunity 
to jamming. In a FH scheme, the jamming power is now spread over a bandwidth of 
B Hz. The signal bandwidth is approximately equal to b Hz. Although the entire signal 
power is present after dehopping, the amount of jamming power present after dehopping 
Is J/b/B). Thus, in a FH system, 


"uy NG) (1.4) 


The processing gain is thus given by B'b. 

In many military applications we would like to deny an adversary the knowledge 
that we are transmitting a signal. That is to say, besides making it difficult for an ad- 
versary to “break” our code once he intercepts our signal, we would like to make it dif- 
ficult for him to detect that there exists a transmitted signal at all. (Obviously a signal 
that is not detected can not be exploited.) To achieve a low probability of detection in 
a communications scheme, the signal power level must be reduced below the noise power 
level. In other words, we need to “bury” the signal in the noise (or hide it behind another 
more powerful already-present signal) so that the adversary will “see” only noise. Such 
a scheme is termed a Low Probability of Detection (LPD) scheme. 

A major disadvantage of the frequency hopping approach lies in the fact that the 
signal can not be reduced below the noise level; in fact, the signal must be maintained 
significantly above the noise level. As a result, frequency hopping schemes are used to 
provide anti-jamming, or low probability of signal detection. 

In some communication environments, multiple transmission paths may arise be- 
tween the transmitter and receiver. This phenomenon, termed multipath, often degrades 
communications since the signals arriving at the receiver via secondary paths 
destructively interfere with the signal on the primary link. In a FH receiver only the 
energy at a particular frequency is passed through the first stage filter at any interval of 


time. If the multipath signals arrive at the receiver with a sufficient time delay, then the 


receiver may have “gone on to the next hop” by that time. In such instances the multi- 
path signals are effectively ignored. This is called the multipath immunity of frequency 
hopping. 

Finally, the utility of the frequency hopping method is limited by the speed with 
which the electrical components can alter frequency without generating excessive noise. 
As of 1990, frequency oscillators were able to vary frequencies at rates of 1 million hops 
per second over a 20 MHz bandwidth [Ref. 3]. 


C. DIRECT SPREADING (DS) 
1. General 

In a direct spreading scheme, the bandwidth spreading 1s accomplished on a 
bit-bv-bit basis. We multiply the data signal v,(r) by a wideband voltage signal v,(r). 
Since multiplication in the time domain corresponds to convolution in the frequency 
domain, the resulting signal’s bandwidth 1s spread. We can then use a matched filter or 
correlator in the receiver to detect the transmitted signal. 

If a matched filter is used in the receiver, an expression can be obtained for the 
system's processing gain. Let E represent the energy per bit in the input signal to the 
receiver, let T represent the bit duration, let HJ be the bandwidth of the spread signal and 
let .V,/2 be the two-sided power spectral density of the noise. Then the input signal 


power 1s given by 








Е 
The input noise power is given by 
y № r 7 
№; = 5 (2H) 2 N,H (1.6) 
and the input SNR is 
Si E ] 
svn, ВЕ (ср) ал 
IN; № TR 


A result from matched filter theory states that the output signal to noise ratio 1s given 
by 


SNR, = E (1.8) 


Combining (1.7) and (1.8) yields 
SAR = SNR OTI) (1.9) 


Thus 2TW is the processing gain for direct spreading using a matched filter. To in- 
crease the processing gain we usually increase IV. Note that there is no theoretical limit 
to the degree of improvement achievable. 

It is a classic myth to believe that direct spreading improves the output SNR in 
receivers where thermal noise limits operation. This is not the case. If the receiver's 
performance 1s limited by thermal noise, then employing spread spectrum will not help. 
Looking at (1.9), we know that utilizing spread spectrum increases WF and hence, we may 
be tempted to conclude, the output signal to noise ratio, SNR, , must increase. How- 
ever, increasing IJ increases the noise bandwidth of the receiver, and hence increases the 


input noise power since 
М = КОИ (1.10) 


where & 1s Boltzmann's constant and T is the effective noise temperature. The increase 
in .V, causes a decrease in SNR The decrease in SNR, counterbalances the increase in 
the /271V) term. 

In jamming limited receivers М, = J —constant. In this case we do get improve- 
ment by employing direct spreading. 

Direct spreading affords us the ability to manage the power in the transmitted 
signal. We can distribute the power over any frequency interval without penalty. We 
can distribute the power over a wide frequency band so that the signal power density 
becomes less than the noise power density. In this way we can “bury” the signal in the 
noise and achieve covert communications. 

To summarize, in the direct spreading technique, we represent each bit of dura- 
tion 7 with a wideband signal v,(r) of duration Г апа bandwidth W. 

2. The Wideband Signal 

The value of the direct spreading scheme depends cnitically on the choice of 
v(t), the wideband signal. The signal must have a good autocorrelation function (ACF) 
since we want a good matched filter output. A "good" ACF is one with a narrow , high 
magnitude main lobe and small side lobes. Such an ACF allows the use of a follow-on 


threshold device with minimum decision error. Additionally, v,(r) should have a 


buildable matched filter and should have a near-white (flat) spectrum to allow for covert 
communications. 
3. Analog Direct Spreading 

In analog direct spreading, the wideband signal v,(r) is analog. The beauty of 
analog direct spreading is that the communication signal can be directly transmitted; 
there 1s no need for a separate modulating carrier. Perhaps the most commonly used 
analog spreading signal is the CHIRP signal. In this case v,(r) varies linearly in fre- 
quency from some initial frequency, f, to some final frequency, fj, over the interval T. 


The formula for the spreading signal 1s 


v1) — A “ЕС + (Ag Loh ја ) (run (1.11) 


To send the bit “1” we may start at f, and terminate at f, (called UP-CHIRP 1f f, 2 f)). 
To send the bit "0" we start at f; and terminate at f; (called DOWN-CHIRP 1f f, 2 f). 
Using Fresnel Integrals it can be shown that the bandwidth of the spread signal is ap- 
proximatelv equal to f; — f. Additionally, the spectrum is approximately constant in the 
interval (f, f) and zero outside this interval. 

The ACF of the analog CHIRP signal is quite good. Figure | on page 7 shows 
one period of the autocorrelation function of v(t) given by (1.11) for the case of 
2TW = 100, T= Isec and f = 50 Hz. The matched filter for CHIRPED signals 1s the 
surface acoustic wave device (SAW) [Ref. 4]. 

4. Digital Spread Spectrum 

In digital direct spreading, the wideband signal v,(t) is digital. Specifically, for 
v(t) we choose a binary sequence. Instead of transmitting a data bit (say a “1”) we 
transmit a binary sequence (say, “1100011”. Each bit in the binary sequence that 
comprises v,(t) is called a chip. As of 1990, the maximum achievable chip rate was 50 
million chips per second [Ref. 3]. 

The digital spreading signal is created by using the output of an z-stage shift 
register with feedback connections. An n-stage shift register can take on, at most, 
2^ —] different states (depending on the feedback connections), so the maximum length 
sequence attainable from the shift register is also equal to 27 21. These maximum length 
sequences, called m-segquences, have the best ACF attainable for sequences of their 
length. The theory of feedback shift register sequences is a well-developed subject, and 
the interested reader is referred to [Ref. 5: Chapter 7]. 
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Figure 1. Autocorrelation Function of CHIRP 


The ACFs of all m-sequences have the same form. The peak value is propor- 
tional to 2^ —1. The sidelobe value is proportional to -1. A normalized m-sequence ACF 
is shown in Figure 2 on page 8. 

Sending data using a spreading sequence (e.g., an m-sequence) 1s called Direct 
Sequence Spread Spectrum (DSSS). To send a "1" we can send one period of the m- 
sequence. To send a ”-1” we can send one period of the inverted m-sequence. For the 
remainder of this thesis we limit our consideration to DSSS systems. 

5. Matched Filters vs. Correlators 
Matched filters for binary sequences usually consist of some variant of a tapped 


delay line (realized as a tapped shift register). The received signal 1s clocked into a shift 





Figure 2. Autocorrelation Function of an M-Sequence 


Source: Couch, L.W. II, Digital and Analog Communication Systems , Macmillan, 1990 


register (of length equal to one m-sequence period) in the receiver. The contents of each 
individual register are tapped, inverted as appropriate and summed to form a discrete 
approximation to the ACF. 

Unfortunately digital matched filters have many deficiencies. A large number 
of taps are required for desired values of processing gain, making the hardware imple- 
mentation (on an integrated circuit (IC) chip) difficult if not impossible. Additionally, 
the chips are “analog” to the extent that they are buried in noise, and so digital hardware 
does not readily lend itself. However we can sull use the ACF property of m-sequences 
and replace the matched filter with a three-port correlator. Unfortunately, but 
unavoidably, data recovery with a correlator requires a receiver-based locally generated 
matching m-sequence, phase synchronized with the received signal’s m-sequence. Of 
course generating the local m-sequence is of no difficulty, but efficiently phase-locking 


it with the received signal has been, and continues to be, a major research concern. 


П. NONCOHERENT SEQUENTIAL ACQUISITION OF M-SEQUENCES 


In this chapter we consider schemes used to acquire m-sequences in direct sequence 


spread spectrum (DSSS) communication systems. 


A. INTRODUCTION 

In order to recover the data from a DSSS communications signal, the signal must 
be despread, 1.e., the data must be extracted by removing the m-sequence. Because 
matched filters are impractical in this application, the autocorrelation property of m- 
sequences must be exploited. One way to remove the m-sequence (and preserve the 
data) is to multiply the incoming signal by a locally generated identical m-sequence 
which is phase synchronized with the incoming signal's m-sequence. Thus, a primary 
function of the receiver is determination of the phase of the m-sequence in the incoming 
signal. The phase is determined by correlating the incoming signal with the m-sequence 
generated in the receiver. The local sequence's phase 1s adjusted until synchronization 
is achieved. Aligning the two m-sequences is accomplished in two stages. The first 
stage, called acquisition, brings the m-sequences into coarse alignment. The following 
stage, called tracking, brings the sequences into precise alignment. 

In this thesis we consider a noncoherent acquisition scheme. (Noncoherent, in this 
case, refers to the modulating carrier, not the m-sequence.) Although coherent schemes 
are simpler and more widely studied, they are not practical because m-sequence acqul- 
sition is usually performed prior to recovering the modulating carrier. 

To acquire the incoming signal, we must examine, using some methodology, the 
various phases that the m-sequence can assume. One methodology would be to simul- 
taneously examine all possible m-sequence positions by utilizing, in the receiver, a sep- 
arate correlator for each potential phase. Such parallel schemes are characterized by 
short acquisition times but complex hardware requirements. On the other hand, we can 
examine each possible m-sequence phase one by one to determine if the local sequence 


is aligned with the incoming sequence. Such serial schemes are characterized by long 


1 This chapter, which presents the model and framework to be used in the subsequent chap- 
ters, is adapted from [Ref. 6 ], a paper that served as the original motivation for this research. This 
chapter, at various times, uses the words, assumptions, notations, equations, methods of organiza- 
tion and conclusions onginally propounded in [Ref. 6]. This matenal is adapted and used with the 
permission of the authors. 


acquisition times but simplified hardware requirements. In this thesis we consider a se- 


rial acquisition scheme. 


В. NONCOHERENT ACQUISITION SCHEMES 

A block diagram of the receiver's acquisition system is shown in Figure 3 on page 
11. The receiver's input signal 1s comprised of the transmitted signal and corrupting 
additive white Gaussian noise (AWGN) with two-sided power spectral density N,/2. The 


input signal at the receiver, assuming no modulation or doppler shifting, is 
(г) + Ара + АТ.) со ог + 0) + м(1) (2.1) 


Where A, is the signal amplitude, a/rj is the m-sequence signal waveform with phase 
iAT, (i taken to be an integer without loss of generality), 7, is the chip interval, A is the 
amount that the phase of the local m-sequence is altered during the acquisition process, 
w, and @ are the frequency and phase of the carrier, and n/t) is the AWGN. The locally 
generated m-sequence waveform 1s described as a(t + ( + y)AT,) where j is an integer and 
I»|« 0.5. The incoming communication signal is multiplied by the locally generated 
signal and then noncoherentlv demodulated. The result Y, 1s checked bv the decision 
processor to determine if the local and incoming m-sequences are lined up to within 
AT./2. If they are not aligned, the local m-sequence phase is altered by AT, and the 
process repeats. If they are aligned, then j must equal i and the tracking circuit 1s then 
employed to reduce y to zero. 


Using the notation in Figure 3, and neglecting double frequency terms, we have 


ni; A 
Ain = | "(да г + (t+ yJAT,) cos(w, Nat = =i TeS, COS EEEN (2:2) 
0 
hie дір , 
Xan + | r(t)a(t+ (t+ yJAT,) sin(w tdi = > TS, SING И (2:3) 
0 
where 
nds 
Nin = | п(ї)а(1 + (j + у)АТ,) cos(w,t)dt (2.4) 
0 
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Figure 3. Noncoherent Serial Acquisition System 


nT, 
Do - | п(ї)а(1 + (} + у)А Т) ѕ1п(о, дах (2.5) 
0 
are independent, zero mean Gaussian random variables with variance o? = nT,N./4 and 


where 


ЛТ. 
SM + | a(t + ¡AT Jal + (j+ Ата (2.6) 
С “0 


For a fixed n, X,, and X, 


,» the inphase and quadrature components, are independent 


Gaussian random variables, each with variance са? and with means 
(4,/2)T,S, cos 8 and (A4,/2)T,S, sin @ respectively. Note that X,, and X,, (and similarly 
X, , and X, ,) are not independent for n x m,. The test statistic for determining alignment 


1S 
Ү,= Х2, + Х2, (2.7) 


The random variable Y, has a probability density function that is noncentral Chi-squared 


with two degrees of freedom [Ref. 1 : pp. 25-29], given by 


mp југа. АУ 
Jy Un) = 22 o Unt An) 2 D P zQ (2.8) 
бп On 





where /4,=[(4,'2)7,S, cos 0]? + [(4,2)T.S, sin 0]? = (А:/4)12:52 апд where L(.) is the 
modified Bessel function of order zero. 

Using Y,, the decision statistic, the acquisition scheme determines if the two m- 
sequence phases are aligned within A7,/2, 1.е., /=1= 0 (we let i =0 to simplify the nota- 
tion), or if they differ by one or more chips, 1e., | - y) - i| ATP; 2 T,. In other words, 


the decision processor decides between the following hypotheses: 


and |yl« ^ 


H, (non-alignment) : | + у| > > 


H, (alignment) · ј=:=0 апд |у|< ^ (2:9) 


Note that the m-sequence phases mav be such that our situation falls between H, and 
H,. In this case 2 < |у+у|< E 

The parameter /, takes on different values under H, and H,. We designate /,, and 
2 л as the worst-case values of 7, under H, and H, respectively. Note that > 72 > 0: 


Using these worst-case parameter values, the likelihood ratio for our test becomes 


Jr, On | Hi) Лар = ОКИ СЗС 


Аму) = a ер( л 
fy, On | Ho) 2e, 0/00, 02) ) 


We now describe the three decision schemes that we will use to resolve between the 
two hypotheses. 
1. Fixed-Sample-Size (FSS) Scheme 
In a fixed-sample-size (FSS) decision scheme, the length of integration is fixed 
a priori and a decision is made based on the resulting test statistic, Y, If the integration 


length is fromr=0 to r= MT. е FSS се сап бе сео и 


>т ова 


FSS: Амр) | (2.11) 


< т say Ну 


where z is the threshold. Since the likelihood ratio is a monotonically increasing func- 


tion of the variable y, [Ref. 6 ], the FSS test is equivalent to the following test: 


>т = Ajj(t) say Н, 


Е, (212) 
<т' = Ajy(t) say Ho 


ESS: sud 


where A;!(.) is the inverse function of A (.). 
2. Sequential Probability Ratio Test (SPRT) 

The sequential probability ratio test (SPRT) consists of testing the likelihood 
ratio against two thresholds [от л = 1,2,3... until one of the thresholds is exceeded. The 
SPRT is optimal for independent Guassian statistics. The length of integration increases 
by one chip each time n increases by one. The threshold that is reached first (upper or 


lower) determines the hypothesis that 1s accepted (H, or H,). The test is described by 


> А бау, 
SIS А (у) | = бау Но (219) 


otherwise, continue to next n 





Since Y, and Y,., are not independent random variables we can not write A,(y,) as a 
product of independent random variables. However, we can place bounds on the Type 
I error (false alarm probability) and the Type II error (miss probability) by using Wald’s 
inequalities [Ref. 7]: 

E? р 


c and Ва--- 








As (2.14) 
where o and $ are the resulting false-alarm and miss probabilities, respectively. If the 
average test length is large, the inequalities can be approximated by equalities. To make 


the real-time SPRT implementation more practical, we can rewrite the test as 


> А(и) = AS (A) say H, 
SERI YM | < В(п) =А, (В) say Hy (2.15) 


otherwise, continue 





where the thresholds, which are functions of n, can be pre-computed. 
3. Truncated Sequential Probability Ratio Test (TSPRT) 
When using the SPRT it is possible that, for certain m-sequence phase dispari- 
ties, the decision scheme may take an excessively long time to decide between H, and 
H,. Such an incident is especially likely when the difference between the phases of the 


two m-sequences is such that it corresponds to a state between H, and M, To avoid a 
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very long test we impose an upper bound on the test length by truncating the SPRT 
(hence the name TSPRT) at a predetermined test duration, and conducting an FSS test 
at that point. The test is described bv 


> А $ау Н, 
Т5р | ifnzn, A 4) | = В say Hg 
otherwise, continue 
я >т зау 
ID Лабу) = Аду) Ж (2.16) 
«t say Hy 


The test is truncated at n = n (provided it has not already terminated). Again we use the 


monotonicity of A,(.) and rewrite the test as 


> А(п) = A; (A) say H, 
TSPRT: ifn<n, v, 4 € Bn) A; (B) say Н, 


otherwise, continue 


^ 2A (т) замени 
Ши=и, у; А (2.17) 
A y 41g 

A key parameter in determining the utility of an acquisition scheme is the time 

that the test takes to acquire the incoming m-sequence. Generally, for a given false- 
alarm and miss probability, the test with the smallest sample size (7) will reach acquisi- 
tion the soonest. The initial phase disparity between the incoming and local 
m-sequences can be regarded as a uniform random variable distributed over the sequence 
period. The hypothesis H, is, therefore, the most probable and so the expected test 
length under H, contributes most significantly to the value of the test acquisition time. 
Hence, it is best to minimize the value of the expected sample size under H, while keep- 


ing the expected sample size to within reasonable values when Hj, does not apply. 


C. DESIGN OF DECISION PARAMETERS 

In this section we determine the values of the thresholds to be used in the three ac- 
quisition schemes. 

Denote the chips of the m-sequence by c,, A =...,0,1,2,..., where ¢, = + Г апа 
the period be N — 2^ — 1, so that c,;y— c, for all k. The m-sequence spro ciii 


waveform can then be written as 
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бо 


а) = У срт: АТ) (2.18) 


К---со 


where p; (1) 1s a rectangular pulse of width T, and height 1, i.e., Рт (1) = 11[070<!< 7, 
and it is zero elsewhere. Letting /=0, we obtain expressions for S, by combining (6) and 
(18): 


n=] 


> cll — ly | Adee + (17 A)r] under H, 


k=0 
5, = < 


Уса о де, Edd under Hy, 
k=0 
п—1 
= #1-1714+ (1314) 2 ска = Sra Under H, 
= | = | (2.19) 
= -8)? cui p 5У 4 = Sao under fly 


k=0 k=0 


where sgn(x) 1s one for x z 0 and -1 for x <0, / is the largest integer no larger than 
С + УЈА апа дӧ = (/ + у)А — 1! Хоте а 1 # 0 апда0<8 <1. 

Defining the per-chip signal-to-noise ratio as 
АШ 


i (2.20) 


Е XE 


we have 





a. 
m | — (SNR)S;, under H, TM 
202 -L(SNR)SÀ, under H, 


The acquisition system does not know the exact values of S,, and S,, in advance; 
however, we can use nominal worst-case values for these parameters when designing the 
thresholds. Simulation results on the schemes suggest that modeling бИ, ~ 0 un- 
der H, and modeling STN and De Jn under H, achieve the desired decision 

k=0 k=0 


errors. Using these approximations, we have 











| z n(SNRYI - | v| А)? = ; under Hi 
я Ра 
m (2,22) 
2c - ^ n0 
d S SMR = E under HH, 
O 


n 


The cumulative distribution function (cdf) of Y, can be written in terms of the Q 


function as 


0 =, 
Ғұ(у,) = . Em 2.23 
nos) ің ү, 5) = 1 ОЦГА а ума ), у, 2 0 E 


where the Q function 1s defined as [Ref. 8] 


OC = | xe € * 5)? p (tx)dx (2.24) 


An iterative algorithm for calculating the Q function is given in [Ref 9 ]. 

We now illustrate design techniques for determining a good choice of thresholds to 
use in the FSS test, the SPRT and the TSPRT such that our decision errors are less than 
our desired limits. We denote « and f as the desired false-alarm and miss probabilities 
respectively. 


First consider the FSS test. Using the Q function, we can write the false alarm and 


Ра = О] мом tom) (2.25) 
Риз = 1 = ON T (2.26) 


From these equations, the values of t’ and M can be obtained by iteratively solving the 


miss probabilities as 


equations simultaneously, such that P, x « and P,,, x f. 

For the SPRT we can use (2.14) to calculate the thresholds by inserting equalities 
in place of the inequalities. 

To devise thresholds for the TSPRT, we split this test into two parts: a SPRT with 


thresholds 4 and B and errors a,,, and f,,,,, and a FSS test with sample size 7i, threshold 
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t and errors o,, and fpe It can be shown (Ref. 10] that the errors of the TSPRT are 
bounded by 


ІЛ 


o трг! C sprt "а J55 


lA 


Върни Ди сй | (227) 


50, to ensure a,,,, < «and f,,,, <= fj, we split these errors into the sums of the errors 
Encso'the SPRT and the errors due to the FSS test. Specifically, we let 


A = б, + бл; = рух + (1 — po)o 
p= В;рн т D = pif uc») (2.28) 


where p, and p, are constants in [0,1]. By setting a,,,, = pya and f,,,, = pip, we can design 
the thresholds 4 and B according to (2.14). Similarly, by setting à, — (1 — p,ja and 
fj, — (1 — p)f, we can design Tand using (2.25) and (2.26). Note that if we set 
men -—0,the ISPRT becomes the FSS test. If we set p, — p, 2 1, the TSPRT becomes 
ШӘОРЕТ. Ғотг values of p, and p, in (0,1) the TSPRT can be considered a mixture of 
peor! and an FSS test. 


D. PERFORMANCE COMPARISONS 

We now present simulation results showing the expected sample size and power 
function for each of our three decision schemes. The expected sample size, called the 
average sample number (ASN), 1s the average number of chips that the test runs before 
deciding between H, and H,. (Note that for the FSS test the ASN is predetermined and 
equal to M). The power function is the probability of accepting H,. For all simulations 
in this thesis, the updating of the local m-sequence generator is half a chip (.e., 
А = 14), the value of y used is 0.5, the m-sequence used has a length of 1023 with prim- 
itive polynomial 1 -- x^ 4- o5 -- x? -- x'?^, and the values of p, and p, (the mixture weights 
for the TSPRT) are both 0.5. 

Figure 4 on page 18 depicts the ASN for the FSS test, SPRT and TSPRT with 
a=f=0.01. Figure 5 on page 18 depicts the power functions for the same three tests. 
All functions are plotted in terms of |j -- y |. When |j - y | 2 0.5 we have hypothesis 
H, (sincej 2 0and y 20.5). When |j - »y| 2 2.0 we have hypothesis H,. The sample 
size M of the FSS test is 258. For all tests the resulting errors are equal to or slightly 


less than our design values chosen for « and fj, so our threshold modeling (2.22) works 


Figure 4. 


Figure 5. 
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rather well. Under H, and H, the SPRT has the smallest ASN. Both the SPRT and the 
TSPRT perform much better than the FSS test under both H, and H,. It has been found 
that the ASN of the SPRT is larger than that for the TSPRT for some range of |j + у | 
between 0.5 and 2.0. (It was found that the ASN of the SPRT can even be larger than 
the sample size of the FSS test when the errors « and ff are chosen to be very small). 
Among all uncertainty phases to be tested in an m-sequence period, the synchroni- 
zation condition occurs only once, the сопішоп 0.5 < |/ + у | < 2.0 occurs at one or a 
few uncertainty phases (depending on the value of A), and the non-synchronization 
condition (#7,) occurs at the remaining phases. So, to reduce the time taken to acquire 
the incoming m-sequence, we should minimize the ASN under /7,. One way to minimize 
the ASN under H, is to utilize the SPRT acquisition scheme, however this results in high 
ASNs for the cases 0.5 «|j -y| « 2.0. It appears that the fastest acquisition will occur 
up cuse the TSPRT with a large value of p, (2.28). This makes the ASN of the TSPRT 
approximately equal to that of the SPRT under H, and results in smaller values of ASN 
under the cases 0.5<|j+y|<2.0. Figure 6 on page 20 shows this comparison for the 


ASN function. Figure 7 on page 20 shows this comparison for the power function. 
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Il]. PERFORMANCE IN THE FADING CHANNEL 


A. GENERAL DESCRIPTION OF FADING 

In a communication system, the received signal may suffer time varving power 
fluctuations. Since these power fluctuations most often represent an attenuation in the 
received signal (as opposed to an amplification), the phenomenon is termed fading. The 
power instability can be considered a result of instantaneous amplification or atten- 
uation of the transmitted signal's amplitude. 

In our previous discussion, the parameters and characteristics of the communi- 
cations channel were considered to be fully known and time-invariant. In reality, espe- 
cially in cases where part of the communications path consists of an unguided medium, 
the channel characteristics vary randomly with time. Utilizing an unguided medium 
(e.g., the atmosphere or ionosphere) as part of the channel unavoidably subjects the 
communications system to the irregular variations which often occur in nature. Channel 
parameters can be altered due to rain and humidity (which absorb microwave energy), 
atmospheric density variations (which refract and reflect electromagnetic waves), 
mountains and other physical obstacles (which reflect or dissipate microwaves), changes 
in the ionosphere’s electron concentration distribution (which affects refraction of High 
Frequency waves), as well as other factors. Semiperiodic and totally random variations 
in the propagation attributes of the channel can be viewed as altering the channel's 
transfer function. 

Although fading arises from a variety of causes, fading phenomena can often be 
modeled as causing multipath distortion. That is, many fading conditions can be mod- 
eled as causing several alternate transmission routes to arise between the transmitter and 
receiver. In addition to the transmission paths which the engineer incorporated into his 
design, additional unplanned transmission paths are excited randomly with time. A sig- 
nal component arriving at the receiver via one path may be out of phase with the signal 
component arriving via a second (and additional) path(s). The arriving signal compo- 
nents interfere (constructively and destructively), resulting in fading. 

Fading channels are classified as flat fading channels and bandwidth-selective chan- 
nels. In a flat fading channel all of the frequencies present in the transmitted signal fade 
in exactly the same manner and the communication signal can be regarded as undis- 


torted in relative shape. The received signal in such cases can be viewed as multiplied 


2] 


by a random variable which accounts for the fading. One example of this type of fading 
occurs in space links passing through a turbulent atmosphere [Ref. 11 : pp.131-132]. 
Flat fading is further categorized as slow fading and fast fading. In slow fading the 
channel variations (which give rise to the signal power fluctuations) are slower than the 
longest period component in the waveform. In such cases we consider the received SIg- 
nal to be multiplied by a constant (although random) amplitude coefficient. In fast 
fading, the channel variations are comparable to the period components in the 
waveform. In such cases we consider the received signal to be multiplied by a fluctuating 
time function which effectively modulates the received signal. 

In bandwidth-selective channels some of the transmitted signal's frequency compo- 
nents are affected to a greater degree than others. Such channel effects are equivalent 
to the insertion of a bandpass filter into the signal transmission path [Ref. 11 : p.132] 
which causes the signal's edge frequencies to have a larger or smaller fade than the center 
frequencies. In such cases a sudden attenuation of just the high frequency components 
may occur, or, alternately, a loss of only low frequency power may be observed. 

It should be noted at this point that fading 1s a major communications obstacle 
which can cause signal attenuations of several tens of decibels. Figure 8 on page 23 
shows the median duration of fast fades as a function of fade depth for a 4 GHz signal 
over a path length of 30-35 miles. 

In the remainder of this paper we narrow our focus to consider only Лаг, slow fading 
channels. 

Many methods have been devised to counter the effects of fading. The performance 
reduction that characterizes fading is caused by the received signal’s amplitude being 
weakened in comparison to its design value. One technique for minimizing the received 
signal’s amphtude fluctuations is to incorporate automatic gain control (AGC) into the 
receiving amplifier (Ref. 12 : p.112]. This is a negative feedback method from control 
theory in which the receiver amplifier output signal is employed to adjust the amplifier’s 
gain. The circuit is designed so that the output level remains steady notwithstanding 
changes in the input amplitude. The major disadvantage of this technique is that the 
channel additive noise amphtude is also increased along with the signal amplitude. 

A better technique to counter fading consists of dividing the communication’s signal 
power between multiple subchannels which fade independently of each other. The 
probability that fading will be extreme in all subchannels simultaneously is very low, so 
if all subchannel outputs are used to reconstruct the signal, we expect better perform- 


ance than over a single channel. The use of multiple transmission subchannels is called 
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Seconds 





Depth of fade - dB below normal 
Figure 8. Fast Fade Duration 


Source: Pierce, J.R., and Pasner, E.C., Introduction To Communication Science And 


Systems, Plenum Press, 1980 


diversity transmission |Ref. 13 : pp.346-350]. Diversity transmission methods fall under 
three general areas [Ref. 14 : pp.632-634]: 


|. space diversity utilizes several receivers in separate geographical areas. The receiver 
outputs are collected to reconstitute the desired signal. When the signal at one 
receiver fades it is hoped that the signal at a distant location (i.e., at the other re- 
ceivers) will be unaffected. 


2. time diversity employs retransmission of the same signal at spaced time intervals. 
It is hoped that if time-varying fading affects a signal at one time, it will not affect 
the signal at a later time. 


3. frequency diversity employs multiple frequency channels to combat the effects of 
frequency selective fading. It is hoped that if the signal suffers fading in one fre- 
quency band, other frequency bands will still be usable. 


B. MATHEMATICAL DESCRIPTION OF FADING 

Thus far we have described fading phenomena in general qualitative terms. Al- 
though we have asserted that in slow, flat fading channels the received signal 1s effec- 
tively multiplied by a constant (albeit random) amplitude coefficient, we have yet to 
grant mathematical justification to this statement. In order to proceed any further with 
the analysis of our receiver in a fading channel, we must determine a suitable math- 
ematical model to describe fading. Since fading phenomena largely depend on unpre- 
dictable variations in the environment, we can not model fading deterministically. 


(Indeed, if it could be modeled deterministically we would then simply account for it by 
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adding an amplification factor in our receiver design and no further consideration would 
be necessary.) The best that can be achieved 1s a statistical description of the amplitude 
coefficient. 

Communication schemes are typically constructed to achieve a preset minimum 
standard of performance. For example, we may design a system so that the probability 
of false receiver synchronization is 0.0] and the probability of failing to detect synchro- 
nization 1s 0.01. Using a statistical description of fading we can simulate our system in 
a fading channel and determine if our minimum standards of performance are attainable. 
If the random channel attenuation factor degrades our system to the point where we fail 
to satisfy our performance standards we can employ three corrective options. 

Our first option is to attempt to modify the receiver design by suitably changing the 
decision thresholds to meet the performance standards. 

Alternately, a second course of action would be to increase our transmitter power 
to a level where it compensates for a worst-case power loss in the fading channel. In this 
Option We mect our performance standards by allowing adequate “design margins.” 

Our third approach would be to relax our performance standards. This option may 
often be the most attractive if we consider other important design variables. As an il- 
lustration of how the third approach may be the most advantageous, consider again the 
system Where we desire that the probability of false receiver synchronization and prob- 
ability of failing to detect synchronization both be 0.01. Suppose further that a simu- 
lation shows that fading reduces both these specifications to 0.1. In our fading channel, 
with our degraded performance standards, it may be the case that the expected time to 
acquire synchronization is one minute. If we employ option one above and redesign the 
receiver, it may be the case that the performance standards are met, but the time to ac- 
quire synchronization is increased to one hour. In such a scenario we would probably 
sacrifice our performance standards for the sake of acquisition time. 

Consider a multipath channel. In addition to the transmission path which the en- 
gineer incorporated into his design, additional unplanned transmission paths of various 
electrical lengths are excited. Portions of the power transmitted at a given time arrive 
“sequentially” at the receiver over each of the sundry transmission paths. Since we are 
restricting our consideration only to flat fading channels, it suffices to examine the 
channel effects on a single frequency component; all other frequencies will be affected 
in a similar way. 

A particular frequency component is transmitted at a certain time. The received 


signals will arrive at the receiver over the various paths. Since the lengths of the paths 
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differ, the received sine waves will (presumably) be out of phase with each other, and 
when combined they will constructively or destructively interfere. The fading channel 
causes the electrical path lengths to vary continually and randomly and this in turn 
causes the phase differences between the component tones to vary. The superposition 
of all the tones (1.e., the received signal) will have an envelope and phase that will vary 
randomly as a result. So, in a fading channel, the received signal is the superposition 
of a number of random phasors. 

Initially assume that each of the random phasors is of a comparable magnitude. 
Consider each of the individual phasors to be resolved into quadrature components. 
Since the fading channel randomly varies the path lengths, the quadrature components 
of one particular random phasor (at a particular time) will be independent of the 
quadrature components of all of the other random phasors. Additionally, the 
quadrature components of any particular phasor are uncorrelated (since the phasor is 
resolved into sine and cosine components). 

Before finding the total received signal, first find the total signal in each of the 
quadrature directions. We do this by simply adding the particular quadrature compo- 
nents of all the random phasors. In doing this we are adding a number of independent 
random variables (from the same underlying distribution) together, and the sum, by the 
central limit theorem [Ref. 15 : p.210], approaches a Gaussian random variable as the 
number of contributing phasors grows large. Thus, at any given time, the two 
quadrature resultants are independent Gaussian random variables. The envelope of the 
received signal will therefore be Rayleigh distributed [Ref. 16 : pp.101-103]. (Note that 
the phase of the received signal is of no concern since we utilize a noncoherent detection 
scheme in the receiver.) 

The received signal, considering all frequencies, will thus be effectively multiplied by 


a coefficient, y, having the probability density function (pdf) 
у -pho 
АЉ) = е (3.1) 
о 


where c?, determined by the nature of the channel, is the variance of each of the 
Gaussian quadrature random variates. 

Since the above discussion utilized the central limit theorem, the question naturally 
arises: how many alternate transmission paths are required for the Rayleigh model to 


be valid? Surprisingly, it has been noted [Ref. 17 : pp.348-351] that as few as six sine 
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waves with independently varying random phases will combine to give a fluctuating re- 


sultant whose envelope closely approaches Rayleigh statistics. See Figure 9. 


Relative Signal Amplitude, dB 
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Figure 9. Envelope Distributions For n Equal-Amplitude Random-Phase 
Phasors: Solid Lines are Rayleigh Distributions. 
Source: Schwartz,M., Bennett, W.R., Stein,S., Communication Systems and Techniques, 
McGraw Hill, 1966, p349. 


In deriving the Rayleigh channel results, we assumed that each of the random 
phasors had a comparable magnitude. [t is more realistic to assume that one of the 
paths, namely the path that we intend to design for (1.e., the no-fading path), is much 
stronger than the other paths. The power in the strong, intended path is termed the di- 
rect component. The total power in all the weak paths is termed the diffuse component. 
The net received signal will then be a combination of a steady tone (again, considering 
a single frequency since the channel is flar) and Gaussian quadrature random variables. 
The envelope statistics for “specular-plus-Rayleigh” fading, called Ricean fading, are 
меп оу ке SE 


HORREA. ono E (3.2) 


с 


where 
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V = instantaneous value of the fading coefficient 
s? = power in the direct component 
20: = power in the diffuse component 


and with the total power having a normalized value of 1. 


: 52 
Petting r= 24 We find that 
NES, 
Лим) 2 29(1 + eT У p OU РА + )) (3.3) 
Note from (3.2) that when the ratio = approaches zero, the pdf approaches the 


Rayleigh pdf (3.1). When r — оо, же have no fading. Figure 10 shows the Ricean 
probability density function (3.3) for various values of r. Note that for large values of 
r the distribution of the fading coefficient is concentrated around unity. 
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Figure 10. — Ricean distribution for various values of r 
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C. TEST STATISTIC DENSITY FUNCTION 

We now determine the probability density function of the receiver's test statistic in 
the presence of a slow, flat fading channel. Since we assume slow fading, y does not 
vary much within one data bit interval, and thus can be assumed constant during the 
acquisition process. 


The received signal through the fading channel 1s given by 


/ 


/ | r(t) = A palt + iAT,) cos(@t + 0) + n(1) (3.4) 


where y 1s the fading random variable with a Ricean probability density function given 


bv (3.2). Referring to our receiver structure (Figure 1), equations (2.2) and (2.3) become 


Ao = 
Xin = VIS, COSA (3.5) 

> а : Е: ~ 
Xan = E WTS, Sin Oa, (3.6) 


where ,V,, and .V,, are still defined by (2.3) and (2.5). As before, the test statistic for 


determining alignment 15 
г ‚2 :2 
m = n + nos (3.7) 


The conditional probability density function of y, given w is noncentral Chi-squared: 


(у 04) 28: 1578 2 » 
е ОМ 295 ly А РИМ on), у, = (3.8) 





| 
TOAS a 
20, 
where, as with no fading, we have /, = (42/4) 72S? and 0? — nT,N,/A. 
We have readily determined the test statistic’s conditional density function given 
V. We are interested in the density function without conditionality. To determine 


f.) we “integrate out” the y dependance: 


ло) = | лолаи 


2 
n 


oo ЕТ? 
= l тоа +АМ) 2e; | A nV Ул V Te Vs 4 
= > 5 е b a e 4 == V 
а 
0 


а с O 
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ра 
2 2 ии 
| е ие 5 2 | Ve ые 2 д Jel E је (3.9) 
0 


2 2 
20,5 On б 


To simplify notation, we combine terms by letting 











¡ba ; с= (3.10) 


and 


A 2 2 2 
d= СЕ пети (3.11) 





Substituting (3.10) and (3.11) into (3.9), we can rewrite the probability density function 


as 


Гог) = ај ре Удар ера (3.12) 


The Bessel function is defined as a power series expansion. Substituting the definition 


of the Bessel function into (3.12) yields 


E = ag 
EET ет Жа 
j=0 


Since integration and summation are both linear operations, the above formula can be 





rewritten as 


б) - ус 20+ Ю+! 2 (3.13) 
me ire Ак!) |: 


j=0 k=0 


It has been derived [Ref. 18 : p. 271] that 
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со ; я 
| рога Да n. (3.14) 
0 


Dae 
Applying (3.14) to (3.13) results in 


pe A ще k)! 
Дин с 
"ME ) бу II 


j=0 k=0 


which 1s rewritten as 





dd — — ас 
аа МП a^ V] ђе МУ 1 (GA. 
мн а (472216) ло o 
j=0 k=0 


The term (jk)! 2 (jt H)fj- 2) (j 3)...(j - Kk) can be written using Appell's symbol 
[Ref. 19 : p. 7] as (¡+ 1,k). Similarly, using this notation, we can rewrite K/ as (1,kJ. 


Adopting this notation, the density function in (3.15) is written as 


ж 
уеэр? 
/=0 


oo 
k=0 


(i +1,k) (b7/4c)* 


m т (3.16) 


The second summation is in the form of the confluent hypergeometric function [Ref. 20 
: рр. 298-320]: 





(2) РИ; 1; 62146) (3.17) 
= 4c (Ат) | 
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Substituting (3.17) into (3.16) yields 


- us a^ V ] с ШО! 
SOn) + 2c )( Ас n a | , 1 ‚Ё /4c) (38) 
J=0 


To further simplify, we express ,F, in terms of a Laguerre Polynomial [Ref. 21 : p. 430] 


as 





ЧОО ee) аа 2. (3.19) 
Emm s 


Substituting (3.19) into (3.18), we find 


d ccu sh ce aa р 2 4 
EA о”) е 


j=0 


where 14.) + L(.) is the Laguerre Polynomial of degree j. Using an identity involving 


Laguerre Polynomials [Ref. 21 : p. 316]: 


k 
) foL) - e AQ) 
К=0 


we can rewrite (3.20) as 


| Яз е И: E 
So) = Ee s ҮЛ ее (3.21) 
: 4c 


where J,(.) is the Bessel function of the first kind of order zero. Using ¿=./— | we have 





2 д : 
Ру») = ze p е“ А (е) (3:22) 


But J,(ix) = h(x), and so we rewrite the probability density function f (y,) as 
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d b* 4c а'4 ab 
fn) = Led er CRY (3.23) 
Substituting the values of a, b, c and d from (3.10) and (3.11) back into (3.23) we have 


our final formula for f (y): 


Е r\(v,/207) + 2,262 
fi.) = ltr Я и 


20: [(1 + г) + 4,/207] (1+7) + 4,20, 


2./r(1 + rl, /207 (2,1207) 
• 1, — (3.24) 
IM 


[his interesting result states that the probability densitv function of y, 1s still noncentral 
Chi-squared. The density function 1s the same form as (2.8) but the parameters 2, апа 


o: are replaced by new vales: 





к IE (3.25a) 
> 5 | а ла бен 
On 4? Оп РР (3.258) 


The cumulative distribution function corresponding to (3.24) 1s 


(2,62) (1-- ғу, 02) 


шала (3.26) 
0702 ЕАС 


A ке 

D. SIMULATION TECHNIQUES 
Now that we have a precise description of the receiver test statistic’s density func- 
tion, We wish to evaluate the receiver’s performance in a fading channel. Three receiver 
characteristics are of primary interest for each acquisition scheme: the actual probability 
of false alarm (PFA), the actual probability of detection (PD) and the number of chips 
needed to make a decision (called the average sample number (ASN) ).2 A fourth char- 


acteristic, acquisition time, which is the average amount of time needed to correctly 


2 Note that the actual probability of false alarm (PFA) and actual probability of detection 
(PD) are distinguished from the desired probability of false alarm (a) and desired probability of 
detection (1 — 8). 
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phase synchronize the incoming and locally generated m-sequences, depends on the 
PFA, PD and ASN. 

There are two general approaches that one may use to determine the PFA, PD and 
ASN for a given fading channel. The first approach is mathematical analysis. This 
method entails determining an equation or formula which allows one to directly compute 
the quantity of interest in terms of other known quantities. Needless to say, closed form 
analytical solutions for the receiver under consideration have eluded, and continue to 
elude, researchers. (The problem, in fact, has not even been solved for the mathemat- 
ically simpler case where the correlators are reset to achieve independent identically dis- 
tributed samples.) 

The second approach used to determine the PFA, PD and ASN for given fading 
conditions 1s simulation. With this method, the threshold values for the various tests are 
stored in computer memory, and a computer program 1s used to simulate the receiver 
operation and the received signal. When the simulation shows that a threshold is ex- 
ceeded, we stop the run and measure the ASN. By running the tests a large number of 
times under various conditions, We can confidently determine the PFA and PD. All re- 
sults presented in this thesis for ASN are based on simulation. 

To illustrate the simulation techniques, we will focus on attempting to determine the 
average ASN (for a particular scheme) with fading. The techniques for determining the 
average PFA and average PD are identical. 

We are interested in determining the effect that fading has on the detection scheme's 
ASN. The average sample number depends on many factors (e.g., noise power spectral 
density, SNR, etc). To isolate the fading effect, we hold all factors stationary except for 
the fading magnitude. So, in other words, we are interested in determining the effects 
of fading on ASN with all other factors remaining constant. 


With the above restriction, the ASN is a function of the fading random variable y, 


ASN = #(4) (3221) 


The expected value of ASN in the presence of fading is thus 


E[ASN] = Elev) = | “aw sway (3.28) 


where f(y) is the Ricean probability density function describing the fading, given by 
(3.2) or (3.3). 
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This is precisely the point where the analysis stops. If we could solve the integral 
in (3.28) we would have the answer we seek. However, since the ASN function (3.27) 
is presently unknown, the integral is not solvable and we must resort to simulation. The 
only purpose in performing the following simulation is to solve (3.28) for E[ASN] (and 
to solve similar equations for E[PD] and E[PFA ). 

1. Monte Carlo Simulation 

We wish to determine a value for the integral in (3.28). The strong law of large 

numbers [Ref. 22 : pp. 88-90] states that for “large enough” m, we can determine 
E[g(v)] (and hence (3.28)) from a sample mean: 


T s) > Elw) (3.29) 
=] 


So, our algorithm for computing Јев) ау is very simple [Ref. 22 : pp. 181-194]: 


l. Generate n independent samples of a Ricean random variable y whose density 
function describes the fading. 


2. Determine g(y,), the ASN for a particular y,, either directly or by simulation for 
Е 


„ә 


The desired result is - за gely). 
=! 


Step 1. above, generating independent samples of a Ricean random variable, is 
of critical importance. If the samples are not Ricean and/or not independent the algo- 
rithm will not be successful. The method used for generating such independent Ricean 
random variables 1s included in the Appendix at the end of this thesis. 

We characterize the severity of fading by specifying a value of r, where r 1s the 
ratio of the power in the direct component to the power in the diffuse component and 
the Ricean density function is described by (3.3). For each Ricean random variate y, 
we run a simulation program to determine the ASN of a test given that value of y. 
We run the simulation program for a large number of independent variates ф,, and 
compute the average ASN by the summation in step 3. above. 

The Monte Carlo method gives us only an approximation of the integral we are 
trying to evaluate. The sample mean that we compute in the third step of the algorithm 
is an estimate of the integral which becomes exact only as n > oo. We are interested in 
determining the mean square error of our estimator. 


Recall that our goal is to evaluate E[ASN] given by 
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СІП #1 (3.30) 


where we designate the exact (albeit unknown) value of this integral by /. Designate the 


estimate of the integral by g,. That is, 


2- 1-2 и) (5-31) 
(=! 


Then the mean of the square of the error of our estimator is given by [Ref. 22 : pp. 
194-196] 


E[(sample mean error)’ ] — E[(g,— ell = Var(g,) (3.32) 


But by using the definition of g,, the variance of g, can be rewritten as 


Var(g,) = var} +) ew) Е E ан) ew} (3.33) 
i=} i=! 


Substituting (3.33) into. (3.32) and using the  shorthand notation 


М = E[(sample mean errory], we have 
У = рага) + (0) + + 800) (3.34) 
, 


Since the variance of the sum of independent random variables is the sum of the vari- 


ances, (3.34) 1s revised: 
У = 1 (матер) + Матео) += + Уаг[е(ф„)]) (3.35) 
H 


But each of the variance terms are equal to each other. That is to say, Var[g(y,)] = 
Var[g(v;)] —..- Var[g(v,)] ^ Var[g(v)] . Using this fact 


У = = (nVarg(/)]) = $ Varie) (3.36) 
n 
But the variance of g(V) is equal to E[g(V)] - (EL[g(v)]? » E[g(v)] - P. Thus 


35 


у = {ЕЕ О ЛЕ (3.37) 


But we know that 


ЕГЕ УР = | g (V) f (U)dy 
0 


so we substitute this quantity into (3.37) and arrive at 


| gSde – 1 
y= 2 (3.38) 


H 


Now, neither of the quantities in the numerator of (3.38) 1s known since we do not have 
an expression for g(w). Thus, we again resort to the law of large numbers to approxi- 


mate each of the numerator terms [Ref. 22 : pp. 196-197]: 


| Pwr Б ТУ (3.39a) 
0 (ее) 


|o "2 #04) (3.395) 


i=] 


So, our final equation for the mean square error for our Monte Carlo estimator is given 
by 


у = {+ : Ye @ -|+ > ava] | (3.40) 


іші 


What is so nice about the Monte Carlo method is that as we compute our estimator we 


can simultaneously compute just how good or bad our estimator is by using (3.40). 
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Again, the exact value of the integral that we are trying to determine is desig- 
nated as J. We would like to be able to place an upper and lower bound on / with a 
certain degree of confidence. Mathematically, we would like to be able to specify a 


number d and then compute two constants c, and c, such that 
В а (3.41) 
Chebyshev's Inequality [Ref. 15: p. 105] states that 


Р[ а-о < X < ptia] > — (3.42) 
I 


where u is the mean of the random variable .Y and ø is its standard deviation. Applying 


the inequalitv to our scenario, we have 


rae uu LEE м (3.43) 


І 


Thus, to choose bounds for (3.41), simply let 


ME (3.444) 





ОЛ К , у (3.445) 


It is known that Chebyshev’s Inequality gives a bound that 1s often pessimistic 
[Ref. 22 : p. 197]. A much tighter bound can be obtained if we use a "large enough” 
number of random variates vy, so that the central limit theorem is applicable. In such 
an event it can be shown [Ref. 22 : p. 198] that 
£l море (3.45) 
JV 





where N(0,1) is the standard normal random variable which we designate as Y. If we 


choose z such that 
d 
I-Fy = PL Y>r] = > (3.46) 


then rearranging (3.43) results in 


37 


р 112,1 <4/У ) = ( | UM 


Ww 





Ш 
= a E 2 = 1- — = [—а (3.47) 
Ку 


As an example of how to apply the above method, suppose we want to deter- 
mine an upper and lower bound for / with 90% certainty. Using (3.41) we would choose 


d = 0.1. Using this value of d, we determine z (from tables} such that 
Е а 
РЕ == p (3.48) 


where Y is the standard normal distribution. Now, with this value of 1, we сап Бе 90% 
certain that / lies between g, — МУ and g, + МУ where g, and v are determined from 
(3.31) and (3.40) respectively. 
2. Numerical Integration 
The Monte Carlo method is but one way to evaluate EL ASN] given by the in- 
tegral in (3.285). A simpler and more direct technique involves estimating the integral 
by using a numerical method. 


The integral we wish to evaluate is given by (3.30) and repeated below: 


ЕСІП (3.49) 
We have a closed form analytical description of f(w), given in (3.24), but we do not 
know the form of g(V). Note that for all v, g(V)f/(v) > 0. 

Our approach is to view the integral as an area. We start by dividing the density 
function of f (V) into a certain number of equal areas. For illustration, we will divide the 
density function into 25 sections, each with an area of 1/25 = 0.04 square units. So the 
first area spans the abscissa from 0 to w,. The second area spans the abscissa from y, 
to w,, and so on. The last area spans the abscissa from фд 10 со. So we can write (3.49) 


as 25 separate integrals: 


205) V V со 
| arto = | елор + [оодар а | о гадар 0.50) 
0 0 Vi Мод 
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Note that the values y, ... V, are easy to determine since f(w) is known. Also note that 
(3.50) 1s still exact. 

Let the midpoints of each of the 25 intervals (0,U,), (V, V3) ... (Y ,,,00) be desig- 
nated in general as y" and specifically as V4, Vs, Vc, ... Vy respectively.3 Note that the 
majority of the mass of f(y) is concentrated in the region between y — 0 and approxi- 
mately y = 3 (depending on the value of r); see Figure 10. Because of this, we expect 
the 25 points V4, Vs, ... Wy to be clustered close together in the region (0,3). Putting this 
another way, we expect each of the first 24 intervals (0,U ), (V, V.) ... (s) to be very 
inn. 

Now we must make an approximation. Suppose over each of the 25 areas, we 
assume g(V) = g(y") where y” 1s the midpoint of the interval (V, V,.,). Since y undergoes 
a very small variation between y, and y,.,, we assume g(V) also undergoes only a verv 
small variation over the interval, and is, essentiallv, constant with a value of g(y’). With 


this approximation, (3.50) can be rewritten as 


[ems] “rr e eto | Prae etis) года зи 
v; : 


0 Ya 


But, by the way that we chose y, v, ... V4, each integral has an area exactly equal to 
№ hus we have 


Јела = e eti) + ea) H «(йул (3.52) 


Our numerical integration algorithm can be summarized as 


1. Divide the probability density function f(y) into X equal areas. Determine the 
midpoint of each area. 


2. For each of the area midpoints, use a simulation program to determine g(y'). 
x 

RE E[ASN] is approximated by = БОО) 
i=] 


A major disadvantage of the numerical integration technique is that we do not 


have bounds on the estimator (or its error) as we did with the Monte Carlo method. 


3 Note that all of the midpoints are readily computed except for Py. Since the 25th integral 
is improper, technically we should have yy = co. To resolve this, we truncate the tail of the 25th 
integral when it goes below a certain very small value, perhaps | x 10-5, and then use the mean 
value theorem for integrals to find an “effective” value for y y. 
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Most numerical integration schemes divide the abscissa into equal areas, and thus have 
well known error bounds [Ref. 23: pp. 285-312]. In our technique the abscissa is spe- 
cifically nor divided up into equal areas, precisely so that we may solve each of the 25 
integrals in (3.51). 

Note that the numerical integration technique is not too unlike the Monte Carlo 
method. For our numerical integration scheme, the points are biased around the mean 
of the density function, as we would expect them to be if they were chosen randomly. 
More precisely, because we divided the density function into a certain number of equal 
areas, each of the y" in the numerical integration scheme are equally likely, and thus we 
can well imagine that thev were generated by our random number generator in the 


Monte Carlo method. 


E. RECEIVER PERFORMANCE WITH FADING 

We now use the Monte Carlo and numerical integration algorithms to examine the 
effects of fading on the receiver acquisition schemes. We begin by examining the per- 
formance of the TSPRT scheme in a fading channel characterized bv r2 ]. Note that 
this is rather severe fading; there is just as much diffuse signal as direct signal. For our 
simulation we use p = p, = 0.5 (2.28), SNR= -10dB, « 2 f — 0.01, and a sequence of 
length 1023. 

1. Expected Results 

Consider again the per-chip signal to noise ratio defined in (2.20) and repeated 
below: 
407, 


SNK = LION (3.53) 


Rewriting (2.21) we have 








SNR = (3.55) 


À 
= E under H, 
| о 20; 
п 
З a under Ho 
Sao 20, 


We know that fading changes the quantities 4, and ø} in accordance with (3.25). 








Therefore, under both 77, and H, it is the case that 
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Е и r 
SNReaping = SA —Ü БЕТТ | (3.56) 
1 +r ++ И 2б; 


In other words, fading has the effect of reducing the SNR bv a factor of 
r/(1-- r4 À,/202). Since 4,f0? 15 alwavs greater than zero, this multiplicative factor will 
always be less than one, and thus the SNR is always effectively reduced below its original 
value, with the exact degree of reduction dependent upon the value of r.4 

Consider again our model used in determining decision thresholds, presented in 


(2.22) and repeated below: 











E 4 
| Z SNR -IIA = ~ under H, 
2” бох 
E (3.57) 
2c? Àno 
n в = > under Hy 
2б, 


We first examine H,. Without fading, our design anticipates /,/202 to be ap- 
proximately equal to the per-chip SNR. With fading, the actual SNR is effectively re- 
duced, so we expect our test to become overdesigned (e.g., overly conservative) under Н,. 
In other words, whereas we accounted for a certain SNR in designing for a low proba- 
bility of false alarm (see (3.57)), the actual SNR is even less than what we designed for, 
thus making our correlator output appear even more "H,-like”. So, we expect fading to 
actually improve PFA (1.е., make it smaller) and lower the ASN under H,. We expect 
our tests, under H,, to perform better as the fading becomes worse. 

Having said this, we quickly add that we expect the effects of fading to be only 
slight under H,. By (3.57), a reduction in effective SNR will mean that our test is over- 
designed by some fixed amount. Since the two m-sequences are not synchronized under 
H,, we should, ideally, detect only noise after the correlation process in the receiver. 
Since we have onlv partial correlation in the decision process, we expect fading to in- 
fluence the test (as described above), but only to a small extent. 

The situation is much different under H, As with H, our design anticipates 
À,/202 to be proportional to the per-chip SNR. However, with fading, our actual effec- 
tive SNR is less than what our test accounts for in attempting to satisfy the desired 


probability of detection. With fading, our signal appears less ^ H,-like". So, we expect 


4 Note that for r —^ co the reduction factor becomes equal to one. This is expected, since 
r= œ represents no fading. 
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fading to cause a deterioration in the receiver's probability of detection, PD, for all ac- 
quisition schemes. 

Although we expect the effects of fading to be slight under H,, we expect the 
effects to be severe under H;. Examining (3.57) we conclude that under H, the test is 
overdesigned by a fixed amount equal to the effective reduction in SNR. Under H,, 
however, the effective SNR reduction is made worse as n increases. Specifically, apply- 
ing (3.56) we see that under H, the SNR is effectively reduced br a constant factom 
(UL +r + SNR]). Under H, the SNR is effectively reduced by the factor 


„с. г ин 
| l+rtn(SNR\1—|y| Ay 


which becomes smaller as n increases. Our scheme under Л, becomes successively more 
underdesigned as n increases. 
2. Monte Carlo Simulation Results 

A 50 run Monte Carlo TSPRT simulation was performed under fading condi- 
tions with r= 1. After the data for all runs were collected, it was arranged in order of 
increasing vy, the fading random variable. Figure 11 on page 43 displays ASN vs. y 
under Н,. Figure 12 on page 43 displays PD vs y under H,. Recall that y = 1 amounts 
to the absence of fading. As y increases bevond у = 1, PD remains high and ASN de- 
creases dramatically. This is expected, since y > I corresponds to a stronger signal (see 
(3.4)) than expected. Our test 1s designed to meet PD=0.99 with no fading; if we have 
“constructive fading” (1.e., > 1) our signal is stronger than what we have designed for 
and therefore fewer chips are needed in the correlation process to meet the 77, threshold. 

As V goes below one, the ASN starts increasing rapidly while PD starts declin- 
ing. The test runs longer in attempting to satisfy PD =0.99, but increasingly decides H, 
by mistake. As у decreases, the H, condition increasingly appears to the receiver to be 
the H, condition. At about wy = 0.55 , PD drops below 0.5 and the decision scheme starts 
deciding H, more often than H,. As w drops below 0.55 the ASN drops as the ISPRT 
decides on H, earlier and earlier. 

Figure 13 on page 44 shows E[ASN] , given by (3.31), vs. the number of runs. 
Note that E[ASN] seems to settle down as the number of runs increases. Under the 
Monte Carlo scheme, E[ASN] would give the exact answer as л — оо. Generally, the 


more runs used, the better the estimator 1s expected to be. 
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Figure 11. TSPRT ASN vs y under H, (Monte Carlo) 
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Figure 12. TSPRT PD vs y under H, (Monte Carlo) 
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Figure 13.  ISPRT E(ASN) vs number of runs: H, 


Figure 14 on page 45 displays ASN vs y under H, Notice that the change in 
ASN over the range of is only about ten chips, a small variation compared to the 
nearly 200 chip fluctuation seen under H, (Figure 11). Figure 15 on page 45 displays 
PFA vs. фр under Hj. Notice that PFA does not go above the design value of « = 0.01 
until y=1.5. This confirms our expectation that fading would have only a small effect 
under Ap 

How confident can we be that our Monte Carlo results are accurate? 
Figure 16 on page 46 shows the mean square error of our ASN estimate (under /7,) vs. 
the number of runs. We expect that as the number of runs increases, the mean square 
error, V, given by (3.40), will decrease, and we see that this is indeed the case. Ideally, 
V —0 as the number of runs — oo. After 50 runs, our Monte Carlo scheme yields 
E[ASN] = 140 with a mean square error ~ 70. Obviously, our results after 50 runs до 
not instill as much confidence as we would have desired. 

A confidence interval calculation is instructive. Suppose we want to determine 
bounds for / with 90% confidence. Using (3.41), we see we require d=0.1. Using (3.48) 
we find r such that P[Y » 7] 2 0.05 where Y is the standard normal distribution. Con- 
sulting a table vields r= 1.645. With g, ^ 140 and V7 70, (3.43) gives our desired con- 


fidence interval: 


ASN 





Figure 14. TSPRT ASN vs y under H, (Monte Carlo) 
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Figure 15. TSPRT PFA vs y under H, (Monte Carlo) 
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Figure 16. MSE of ASN estimate for ISPRT under Н, 


P[126 < I < 154] > 90% 


In other words, we can be assured that E[ASN] under H, 1s within the 28 chip band 
(126-154) with 90% confidence. 
3. Numerical Integration Simulation Results 

We reran the TSPRT simulation under г = 1 fading conditions by dividing the 
Ricean density into 23 equal areas. Figure 17 on page 47 shows the ASN vs y for H, 
Note that these results, and the results to follow, closely resemble those obtained 
through Monte Carlo simulation. Figure 18 on page 47 displays PD vs. y under Hj. 
Figure 19 on page 48 and Figure 20 on page 48 show ASN and PFA, respectively, vs. 
w under H,. Figure 21 on page 49 includes data for PD vs. y for points between Н, апа 
H 

How confident can we be that our numerical integration results are accurate? 
Recall that there is no precise way to determine confidence intervals for our integration 
scheme. To gain insight into how stable our results are, we reran the integration scheme 
but divided the density function into 50 areas instead of 25. We found that, to the 
nearest chip, the value of E[ASN] under H, for 50 areas was identical to the value ob- 


tained for 25 areas. Similarly, the value of E[PD] using 50 areas was found to be equal 


46 


LEGEND 
о = САММА=О.4 
© = GAMMA=0.5 


ASN 





0 025 050 0.75 1 125 150 1.75 2 
РЫ 


Figure 17. TSPRT ASN vs y under H, (Numerical Integration) 
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Figure 18. TSPRT PD vs y under H, (Numerical Integration) 
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Figure 19. TSPRT ASN vs y under H, (Numerical Integration) 
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Figure 20. TSPRT PFA vs w under H, (Numerical Integration) 
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Figure 21. TSPRT PD vs y (Numerical Integration) 


to the value obtained with 25 areas to within three decimal places. We conclude that 
our numerical integration scheme using 25 areas is accurate and, at this point, we 
abandon the Monte Carlo method for all future simulations in this thesis. 

We summarize our results for the TSPRT with r2 1.. Without fading, our ex- 
pected value for the probability of detection under H, is 0.992. With fading, the expected 
value for the probability of detection dropped to 0.762. Our expected value of ASN in- 
creased from 138 without fading to 156 with fading. 

Under H, our expected value of probability of false alarm decreased from the 
no-fading value of 0.01 to a value of 0.007 with fading. The ASN under H, decreased 
from 158 without fading to 155 with fading. 

Figure 22 on page 50 shows the TSPRT power vs. |j 4- y | with and without 
fading. Recall that |j - y | 2 0.5 corresponds to МН, while |j -- y | 2 2.0 corresponds to 
H, Figure 23 on page 50 summarizes the fading results for ASN vs. |j 4 у |. 

4. Simulation For Various Fading Conditions 

We obtained simulation results for the TSPRT with fading for the additional 
cases of r=0, 10 and 20. In all cases the simulation technique consisted of numerical 
integration with the density function divided into 25 equal areas. Figure 24 on page 51 


summarizes the results for ASN vs. |j - y | while Figure 25 on page 51 summarizes the 
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Figure 22. TSPRT Power vs |j + 


ASN 


LEGEND 
GC = NO FADING 
© = FADING R=10.0 





0.50 0.75 1.50 1.75 4 


1.25 
| +САММА | 


Figure 23. TSPRT ASN vs |j+y| 
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Figure 24. TSPRT ASN vs] + y | (various г) 
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Figure 25. TSPRT PD vs |j 4 y | (various г) 
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results for PD vs. |j - y| Note that the test performs better as r increases, as expected. 
Additionally, it is seen that fading has very little effect on the H, hypothesis, as expected. 

We then simulated the FSS test under identical conditions for r=0,1,10 and 20. 
The results for PD vs. |j - y | are summarized in Figure 26. (Note that for the FSS test, 
the ASN is a constant independent of fading.) 
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Figure 26. FSS PD vs|j+ y| (various r) 


F. THRESHOLD ADJUSTMENT 

We have seen that fading causes no deterioration of the receiver performance under 
H, but generates a considerable degradation under H, Under worst possible fading, the 
Rayleigh case, the probability of detection is reduced from 0.992 to 0.724 for the TSPRT 
and from 0.990 to 0.721 for the FSS test. 

Since we know the exact form of the test statistic's probability density function with 
fading (given by (3.24)), we can redesign our decision processor to account for the fad- 
ing. Designating 4, as 2,, under H,, and similarly designating /, as 4, under A,, we write 


the likelihood ratio as 


| +” + (2,203) 2.0.289 2 11205 
Ann) СД M 
| Фу + (2, 1/26,) | + у + Ла 267 l+r+24, 1/20, 
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Regrouping terms, we can isolate all terms with the quantity у„/2о; into a test statistic 


and compare this test statistic to two thresholds: 
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/202 аге given by (3.57) and A and B are determined by 


п,0! 


the procedure detailed in Section C of Chapter II. For the FSS test A = B and our two 


where, as before, 2, /202 апі 2 


thresholds are equal. 

Note that the new test statistic, as well as the thresholds, depends critically on the 
value of r. Adjusting the thresholds to counter the fading effect requires that we some- 
how measure quantitatively the value of r prior to profiting from the updated design. 
In other words, if we suspect that fading has caused our receiver performance to deteri- 
orate, we must determine the degree of fading, r, prior to adjusting the thresholds. A 


detailed study of techniques that may be used to determine the fade depth is beyond the 
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scope of this thesis, but the following methods, individually or in combination, might 


be considered: 


1. If the fading 1s due to multipath, we can employ a highly directive antenna to scan 
the received signal strength in azimuth and elevation. Using this information, we 
determine how much of the received signal power is arriving from directions outside 
the design path. 


2. If the fading 1s due to ionospheric effects, we can use sounding techniques to 
measure the layer alterations, or we may avail ourselves of the historical informa- 
tion already gathered concerning the ionosphere. Figure 27 on page 55 [Ref. 24 ] 
shows the fade depth in dB for a typical HF channel at mid-latitudes, over a 200km 
path. We can use (3.56) to convert this information into an estimate for the value 
ofr. 


3. Prior to initiating communications, we can conduct an extensive channel survey to 
measure and assemble information on signal performance under varying conditions 
(e.g., rain, Magnetic storms, sun spot activity, channel traffic, etc). If similar con- 
ditions arise again once the channel is in use, we will have a baseline estimate of 
the amount of fading to expect. 

Howsoever determined, we assume for the remainder of this discussion that we are able 
to gauge the fade depth and calculate r. 

Figure 28 on page 56 displays the degree to which the thresholds in (3.59) vary with 
fading for the TSPRT with SNR=0.1, length= 1023 and p, = p, =0.5. Note that the 
“original” thresholds correspond to r > oo, or no fading. As r decreases, fading becomes 
Worse and the test runs for a longer time prior to truncation. This makes intuitive sense 
when one considers that fading effectively reduces the received signal's SNR, and so we 
must integrate more m-sequence chips to accumulate enough signal to make a decision. 
Note that if the fading becomes too severe, the TSPRT truncation point becomes greater 
than the sequence length. 

Figure 29 on page 57 shows three typical test statistic trajectories and their inter- 
action with the thresholds for the case of r 2 10. Note that for the A, sample path with 
y = 0.46 (severe fading) the decision processor runs all the way to test truncation prior 
to resolving between the two hypotheses. Note that as the fading worsens from r = оо, 
100, 20 to 10 the test truncation increase from 306 to 332, 461 and 710 chips respectively. 
The truncation point for r = | (not shown) is an extraordinary length of 15,221 chips. 

We reran the TSPRT with simulated fading, but with the thresholds adjusted in ac- 
cordance with (3.59), for the cases of r=10 and r=20. As before, SNR=0.1, 
« = В —0.01, p, and p, 7 0.5 and the sequence length is 1023, The resultant power vs. 
Jt y, the phase offset, is displayed in Figure 30 on page 58. The power curves with no 


adjustments to the thresholds, presented in Figure 25, are repeated for purposes of 
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Figure 27. Fade Depth For Typical HF Link 
Source: Maslin,N., HF Communications: A Systems Approach, Plenum Press, 1987 


comparison. We see that the TSPRT is fully restored; adjusting the thresholds raises the 
probability of detection to greater than 0.99. Threshold adjustment also retains the 
probability of false alarm below 0.01, as desired. 

While we may take satisfaction in restoring the TSPRT power, there is a severe price 
to be paid in employing our decision processor. Figure 31 on page 59 displays the 
TSPRT average sample number (the number of chips the decision circuit must correlate 
prior to reaching a decision) vs |j -- y |. Note that the ASN has increased drastically. 

Recall that among all uncertainty phases to be tested in an m-sequence period, the 
synchronization condition occurs only once, the condition 0.5 «|j - y | « 2.0 occurs at 
one or a few uncertainty phases (depending on the value of A), and the non- 
synchronization condition (H,) occurs at all the remaining phases. For our TSPRT 
simulation, the condition 0.5 « |j + y | « 2.0 occurs at three uncertainty phases while the 
H, condition occurs at 2042 uncertainty phases. To reduce the amount of time taken to 
acquire the incoming m-sequence it is critical that the ASN be minimized under H, With 
no adjustments to the thresholds, E[ASN] under H, is 158 chips. When the thresholds 
are adjusted, E[ASN] increases to 215 chips for r 2 20 and 308 chips for r=10. In 
summary, the price to pay for meeting the probability of detection is a large increase in 
E[ASN] under M, Similar results were observed for the FSS test: the test length in- 
creased, favoring the probability of detection while sacrificing acquisition time. For 
a =f E 0.01 and SNR 7 -10dB, the required FSS sample size increases from 258 (for no 
fading) to 361, 509, 6617 and 8944 when fading is characterized by r- 20, 10, 1 and O, 


respectively. 


55 


% LEGEND 
o = LOW THRESHOLD, ORIGINAL 
à = HIGH THRESHOLD, ORIGINAL 
+ = LOW THRESHOLD, А=10 
x = HIGH THRESHOLD, R=10 


28 © = LOW THRESHOLD, R=20 
v = HIGH THRESHOLD, R=20 
u = LOW THRESHOLD, R=100 
x m HIGH THRESHOLD, R=100 





20 


18 x 


THRESHOLD 


0 100 200 300 400 500 600 700 800 


Figure 28. TSPRT Thresholds Adjusted For Fading 
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Figure 29. Typical Trajectories For r= 10 
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Figure 30. TSPRT PD Fading Performance 
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Figure 31. FSPRT ASN Fading Performance 
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IV. PERFORMANCE WITH MODULATION 


A. GENERAL DESCRIPTION OF MODULATION 
Up to this point, we have only discussed the acquisition of a communications signal 


that does not carry data. All the prior analysis assumed that our received signal was 
r(t) = A a(t + ГАТ) соѕ(о, + 6) + n(1) (4.1) 


where A, is the signal amplitude, a/t) is the m-sequence signal waveform with phase 
iAT, (i taken to be an integer without loss of generality), J, is the chip interval, A is the 
amount that the phase of the local m-sequence 1s altered during the acquisition process, 
о), апд 6 аге («ће frequency and phase of the carrier, and n/t) is the AWGN. 

In this chapter we examine the performance of our detection schemes in the pres- 
ence of data. By far the most commonly used modulation scheme is binary phase shift 
keving (BPSK). For our discussion we assume that the data rate is equal to the m- 
sequence chip rate divided bv the sequence length. In other words, we assume we send 
one bit of information for each period of the m-sequence. 

In the BPSK modulation scheme, to communicate a “1” we transmit one period of 
our chosen m-sequence. To communicate а “-1” we transmit the complement of the 
m-sequence. The m-sequence can thus be modulated by the data bv using a single 
EXCLUSIVE-NOR gate, with the data bit stream and m-sequence chip stream as in- 
puts. The modulated m-sequence is then multiplied by a high frequency carrier. for 


transmission. Representing the data stream as (t), our received signal is now 
r(t) = Ajd(i + ¡AT Ja(: + ¡AT.) cos(w,t + 0) + n(1) (4.2) 


By examining Figure 3 and (2.6), we see that modulation may present a considerable 
problem for our receiver. Note that in the correlation process we multiply our locally 
generated m-sequence by the incoming m-sequence, and then integrate this result. If 
we transmit the data "1^ during the duration of the correlation process, then (2.6) will 
apply just as before. If we transmit a "- 1" during the duration of the correlation process, 
(2.6) will be preceded by a negative sign, which is subsequently removed by the squaring 
process in the decision statistic (2.7). No matter what data we send (“-1” or “1”), the 
receiver performs exactly as it did before, so long as the data remains unchanged during 


the correlation process. 
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The problem arises when a data boundary is encountered, and the data changes sign 
(“-1” to “1” or “1” to “-1”) during the correlation process. Under the ff, condition the 
correlator output builds up in magnitude until a decision is made. When a data sign 
change is encountered, the correlator output starts building in the opposite direction, 
wiping out the previous gains. As an example, suppose we start by sending a ^-1". In 
the /7, condition, the quantity S, , given in (2.6), will start building “up” in the negative 
direction. (The fact that S, becomes more and more negative 1s no concern since we are 
interested in the square of the quantity $,.) Now, if the data bit changes to "1^ during 
the correlation process, the per-chip S, becomes positive. This positive quantity adds to 
(and thus diminishes) the previously accumulated negative running total. 

It is true that we can utilize our communications scheme in such a way that we first 
acquire the m-sequence and then, after acquisition, start sending data. Such a method 
is, in effect, utilizing a data preamble of “1111...” to assist in m-sequence acquisition. 
We are interested in seeing how our acquisition schemes perform when data is trans- 
mitted from the outset, without a “1111...” preamble. Knowledge of performance with 
modulation is useful for the following reasons: 


1. If we decide to transmit a signal without data for a certain fixed period of time in 
order to assist receiver acquisition, then ıt is likely that in many instances the re- 
ceiver will be synchronized for a certain time before actual data transmission com- 
mences. Suppose, for example, that we transmit our signal without data for one 
minute in order to allow the receiver to acquire the m-sequence. Consider that 
when two m-sequences are synchronized, such a condition can be detected in about 
one second. Suppose further that it just so happens that, by chance, the receiver's 
m-sequence is immediately in phase with the incoming m-sequence at the start of 
the acquisition process, 1.e., suppose that we are immediately in the H, state and 
we acquire in about one second. Then, in such a scenario, we “waste” 59 seconds; 
our receiver is “ready to receive” but no data is sent for a considerable time. If, on 
the other hand, we send data immediately, we expect that it may take longer to 
acquire but, once the receiver is synchronized, data transfer immediately com- 
mences. Now, it is true that if we send data immediately, some data sent before 
acquisition might be lost, and the time needed to acquire the signal might be 
greater. On the other hand, more data might be sent because we do not “waste” 
channel time. Knowledge of the performance with modulation may help us resolve 
this tradeoff. 


2. The spreading waveform a(r) is periodic and therefore is characterized by a line 
spectrum in the frequency domain. Multiplying a(t) by the cosine carrier also re- 
sults in a quantity that has a line spectrum in the frequency domain. If we are 
trying to conceal our communications from a potential adversary, we would prefer 
that our signal had a continuous spectrum instead of a more easily noticeable line 
spectrum. For this reason, in secure DSSS communications, it 15 not advisable to 
send preambles of all ones or all zeros. Once we start sending data, the random 
(1.e., non-periodic) nature of the data results in a continuous spectrum for our 
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communications signal. We conclude that sending data immediately affords us a 
more secure communication signal than one that begins with a “1111...” preamble. 


3. No communications scheme is perfect, and we can expect that our scheme will, at 
some time, lose synchronization after data transfer has begun. Our system may 
be operating for a period of time, after which the locally generated and received 
m-sequences fall out of synchronization and make data recovery impossible. 
Causes for losing lock include drift in the oscillators in the transmitter and/or re- 
ceiver, equipment failure and operator error. After such an incident has occurred, 
we must go through the acquisition process all over again. During such an event 
the transmitter will (most likely) not know that the receiver has lost lock, and so 
the transmitter will not cease sending data. The receiver will have to reacquire in 
the presence of data modulation. 

B. TEST STATISTIC DENSITY FUNCTION 

We now determine the probability density function of the receiver's test statistic in 
the presence of data modulation. We make one important presupposition: we assume 
that in every case the number of chips used by the decision processor is less than the 
m-sequence length. This 1s the same as saving that we assume that we will encounter 
at most one data boundary during the correlation process. 

We assume that the test statistic under H, is unchanged when data modulation is 
introduced. When the sequences are out of phase we have, essentially, only noise at the 
correlator outputs and so we do not expect modulation to have much of an effect. Our 
confidence in this assumption is further reinforced by the fading simulation results which 
showed that fading had virtually no effect on the H, hypothesis. Under H,, with modu- 
lation, the pdf is still given by (2.8) with /,/2o?= SNR. Similarly, the cumulative dis- 
tribution function is still given by (2.23). 

We now consider the H, hypothesis. When we transmit data there will be a data 
modulation boundary every N chips, where N is the number of chips in one period of the 
m-sequence. We assume that the location of the data boundary is uniformly distributed 
over the m-sequence length. Specifically, if we start the receiver correlation process at 
t = Q, then, if i is the number of chips integrated prior to encountering the data bound- 
ary, iis a discrete uniform random variable on {0,1,2... ,N — 1}. 

For each decision scheme, the correlation process will run for a certain (unknown) 
number of chips prior to exceeding a decision threshold. Let the number of chips in the 
integration interval for a particular test be n. In other words, our correlator runs from 
t= 0 and makes a decision at r=nT,. The probability that the modulation boundary 
will nor appear in the integration interval is thus 1 — 1/N. Now, suppose that a modu- 


lation boundary does appear in the integration interval, i.e., at some point in the first 
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n— 1 chips of the correlation process. Let the data prior to the boundary be designated 
d, and the data after the modulation boundary be designated d, where 4,4 є {1, – 1). 
Since the data is random, the probability that a sign change occurs at the data boundary 
is 2 :P(dj d) ^. Alternatelv, the probability that no sign change occurs at the 
modulation boundary is also 2. Combining this with the fact that a modulation 


boundary will only appear in the integration interval with probability n/N, we conclude 


Prob(modulation has no effect) = (1 Б A + E т = = Е (4.3) 

When modulation has no effect, the pdf of the test statistic is, as before, given by 
(2.8) with /,/202 given by n(SNR)(1 — | y | A)’, in accordance with (2.22). 

Now, we examine the one remaining case: a data boundary does occur and, at the 
boundary, the data changes sign. Reviewing Chapter II, we see that the only term in the 
description of the receiver decision processor which has changed is S,, given formerly by 
(2.6), and previously modeled by (2.19). Using the same model as in Chapter II, we 


obtain a new expression for S, in the presence of data modulation: 


і-! 


5, = У – |у| Ајс, + 190145655] 


k=0 
n—1 
+ а) а - 11 Аск + || Асо) (4.4) 
k=i 
Using the fact that d = — d,, we find 
im] п--1 
D. == да a п)(1 аа | у | А) - | y | Аана ai y] 
k=0 k=i 


As before, we assume each of the summations is approximately equal to zero, so it fol- 


lows that 
S, = do2i 11714) (4.5) 


Since (2.21) still applies, we have 
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where /,,/202 1s equal to n(SNR)(1 —| y | A). 
Combining (4.3) and (4.6) with (2.8), we determine the closed form probability den- 


sity function for the test statistic under the 77, hypothesis as 
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The corresponding cumulative distribution function, using the Q function defined in 


(2.24), 1s given by 
Fy On) = (-3ғ/1-9 = Ре) 
n=l - 
dTa a «m 
i=0 


If we choose to work with a normalized random variable z, where z = y,/o? , then (4.7) 




















is rewritten as 
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We see from (4.9) that the value of z will depend on the number of chips, n, in the inte- 


gration interval. The expected value of z for a given value of n is 


ЕГ] = (1 E je + 2n(SNR)(1 — | y | Ay’) 


п—1 


+ == У (2+[ 28 ] опузхвја 1714) (4.10) 
і=0 


The quantity E[z] 1s plotted vs. 1 in Figure 32 on page 66 


C. RECEIVER PERFORMANCE WITH MODULATION 
We now use the numerical integration algorithm to examine the effects of fading on 
the receiver acquisition schemes. For all simulations we assume SNR = -10dB and the 
m-sequence is of period 1023 chips. Additionally, for the TSPRT simulations, we as- 
sume p, = р = 0.5. 
1. Expected Results 
We expect that modulation will not have a severe effect on the performance 
(characterized by PFA, PD and ASN) of our receiver acquisition schemes. When a 
modulation boundary is encountered, half the time it will not represent a data sign 
change and thus will be transparent to the receiver. 
To reach a more quantitative estimate of the effects of modulation, we can ex- 
amine the mean value of the quantity 4,/2e2 considering its value when a data change is 


encountered (given by (4.6)) and when no data change is encountered (given by (2.22)): 





dn, ; 
Be» (1-4 )ntsNRa – 1514) 
б, 
п—1 
fut DS _ 2 
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This equation can be further simplified as 
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Figure 32. Expected Value of Z vs n 
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Comparing (4.12) to the original value of 4,,/2e? given by (2.22), we see that 


modulation has the effect of reducing the SNR by approximately 


| n 1 
¢ ЗА А ) 


For large values of NV, the reduction factor is small. In the worst case, n= N and the 
SNR is reduced by 1/3, or 1.8dB. 
2. Numerical Integration Simulation Results: FSS Test 
We examine the effects of modulation on the FSS test assuming the original 
receiver structure. In other words, without accounting for the fact that modulation al- 


ters the test statistic’s density function, we will determine the receiver performance. 
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We first simulate the case where the data boundary is encountered and repres- 
ents a change in the data sign (-1 to +1 or +1 to-1). If no data boundary is encount- 
ered, or if one 1s encountered and it does not represent a sign change, the modulation is 
transparent to the receiver and the resulting PD and PFA are as determined in Chapter 
II. (Recall that the ASN for the FSS test is predetermined by the desired probabilities 
of false alarm and miss: « and f.) 

Let us examine a particular case in detail. Suppose we choose « 2 f/ 20.01. In 
this event, the test termination (found by iteratively solving (2.25) and (2.26)) turns out 
to be 258 chips. We then run the FSS test for each possible value of the data modu- 
lation boundary ; (1.e., we run the test for i = 0,1,2, ... 257). 

Figure 33 on page 68 displavs the actual probability of detection vs. the mod- 
ulation chip boundary under the Н, hypothesis. We see that if the modulation boundary 
is encountered early in the integration process, the test can still reasonably recover and 
yield a relatively high PD (although significantly less than 1 — f). Similarly, if the 
modulation boundary is encountered very late, the integration may have built up suff- 
ciently to result in a relatively high PD. The worst case occurs when the modulation 
boundary is midway in the integration process. In such an event, we see from (4.5) that 
S,=0 and the decision scheme decides on HT. 

Figure 34 on page 68 shows the actual probability of false alarm vs modulation 
chip boundary under the H, hypothesis. We see that modulation does not seem to have 
a major effect on PFA. Figure 35 on page 69 shows PD vs the modulation boundarv 
for some points between H, and A. 

Figure 36 on page 69 summarizes the results for the FSS test power vs |j + у | 
in a scheme employing data modulation. With no data modulation, PD is 0.990. When 
data modulation is present PD drops to 0.921. The PFA remains at the same original 
value, to three decimal places, when modulation 1s present. 

For our given test, with our given sequence length, the modulation effect does 
not seem very severe. As a matter of fact, simulation shows that for this particular set 
of conditions (SNR, a, £,N, etc), modulation causes the same performance degradation 
as fading with r = 10. 

We now evaluate the performance of the FSS test with a variety of values of 
« and f. Specifically, we wish to evaluate the effects of modulation for 
0.01 < œ,ß < 0.20. 

To accurately evaluate the effects of modulation on the FSS test, we must first 


ascertain how the test performs without modulation for all cases of interest; in other 
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words we must determine baseline values of PD, PFA and test termination M for all a 
and f which satisfy 0.01 € 0,65 < 0.2. Figure 37 on page 71 depicts the calculated test 
termination point (assuming no modulation) vs. the design probability of false alarm 
(x) and design probability of miss (f). The termination point is calculated using the 
methods of Chapter II. We note that the termination point seems to be slightly more 
dependent on fj than on a . 

Figure 38 on page 72 displavs the actual probability of miss, 1-PD, hereafter 
designated as PM, vs. a and f in the absence of data modulation. The test performs as 
expected. For a given o, PM is approximately equal to f. For a given fj, РМ 15 ар- 
proximately constant as a varies. Figure 39 on page 73 displays the results for PFA vs. 
c and fj with no modulation. 

We now run the FSS test for all cases of interest (0.01 < a,f < 0.2) in the pres- 
ence of data modulation. The test termination point vs. a and £ is still as shown in 
Figure 37 since , at this point, we have not adjusted the thresholds. Figure 40 on page 
74 shows the PM as a function of the design miss and false alarm probabilities. Com- 
pare this figure to Figure 38. Figure 41 on page 75 shows the variation of PM with « 
for a fixed value of f. The results can be explained by noting that as we relax our design 
probability of false alarm, the FSS terminates sooner (see Figure 37), and thus there is 
less chance that a modulation boundarv will be encountered in a given integration in- 
terval. So, as we increase a, modulation tends to become more transparent to the re- 
ceiver, and PM moves closer to ff. Figure 42 on page 75 shows the variation of PM with 
p for a fixed a. The curves are linear (as desired), but they are offset higher than the 
ideal characteristic, Which would be a line of slope | passing through the origin. Note 
that the degree to which the curves are offset from the ideal curve decreases as a in- 
creases. As we increase a, the test termination point decreases and modulation becomes 
more transparent to the receiver. 

Figure 43 on page 76 displays the actual probability of false alarm vs a and f. 
Note that this curve for PFA is almost identical to the results obtained when no modu- 
lation is present, shown in Figure 39. Modulation seems to have almost a negligible ef- 
fect on performance under H,. 

To summarize our results for the FSS test, we note that modulation worsens the 
actual miss probability (Figure 38) and has a negligible effect on the actual false alarm 
probability. The degree to which the actual miss probability (with modulation) exceeds 
the desired miss probability depends on the design values chosen for « and f. As we 


increase à. and/or ff the test terminates earlier (Figure 37) and therefore modulation has 
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Figure 37. Calculated FSS Termination Point 


less of an effect on receiver performance. As we increase a and/or £, the actual proba- 
bility of miss more closely approaches its design value. 
3. Numerical Integration Results: TSPRT 
We examine the effects of modulation on the TSPRT, assuming the original re- 
ceiver structure. Again, we designate i as the data modulation boundary, and we assume 
iis a uniform discrete random variable on (0,1,2 ... N — 1). 
When we utilized the numerical integration algorithm to simulate modulation 


for the FSS test, we essentially divided the density function of ¡into 1023 equal areas. 
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Figure 38. Actual FSS Probability of Miss 


Since the TSPRT simulation program takes a considerably longer amount of time to 
execute, we can not divide the density function of i into so many equal segments. In- 
stead, for all TSPRT simulations, we divide the uniform pdf into X equal areas, such that 
approximately twenty modulation boundaries may appear at equallv spaced points in 
the integration interval. An example will illustrate our approach. For the TSPRT with 
а = fj — 0.01 , the test is truncated after 306 chips have been integrated (if the test has 


not already ended). We divide the uniform density of i into 64 equal areas of 16 chips 
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Figure 39. Actual FSS Probability of False Alarm 


each.5 We then run the TSPRT with each modulation boundary: 8, 24, 40, 56...296; 19 
runs in all for this case. Our final results for PD, PM ASN and PFA will represent val- 


ues averaged over all the modulation boundaries. 


5 Since the m-sequence length is always an odd number, we either pad the sequence with an 
extra chip, or subtract a chip, before dividing the pdf of i into a certain number of equal areas, each 
with an integer number of chips. 
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Figure 40. Actual FSS Miss Probability With Modulation 


We start by examining in detail an arbitrary case of the TSPRT. We choose a 
— 0.04 and f = 0.14. In such a communications scheme, we do not want to miss the 
H, condition more than once in seven times and we do not want to falsely declare the 
H, condition more than once in 25 opportunities. The truncation point occurs after 162 
chips are integrated. We examine the effects of modulation with no threshold adjust- 
ment. The pdf of i 1s divided up such that there are 21 equally spaced potential modu- 


lation data boundaries within the integration interval. 
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Figure 43. Actual FSS False Alarm Probability With Modulation 


Figure 44 on page 77 displays the average sample number vs. modulation chip 
boundary under H,. Figure 45 on page 77 shows the actual probability of detection vs. 
the modulation chip boundary under H, . We see that if the modulation boundary 15 
encountered very early, the test will run longer and attempt to decide the H, hypothesis. 
If the modulation boundary is encountered after the correlation process is well under- 
way, PD drops dramatically. The ASN also drops since the test simply decides H, at an 
earlier and earlier point. The PD reaches a minimum if the data modulation boundary 


is encountered at approximately the 45th chip into the receiver correlation process. If 
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the modulation boundary 15 encountered after this point, PD and ASN start increasing; 
the later the modulation boundary 15 encountered, the less damage done to the inte- 
gration result that has built up prior to the boundary. 

Figure 45 depicts the variation of PD vs. modulation chip for the TSPRT, while 
Figure 33 shows the same information for the FSS test. The minimum value of PD oc- 
curs in the FSS test at the midpoint between chip 0 and the termination chip. The 
minimum value of PD occurs in the TSPRT at a point far to the left of the midpoint 
between chip 0 and the test truncation value. This difference in the two curves arises 
from the fact that, whereas the FSS always runs the whole span from chip 0 to the ter- 
mination chip, the TSPRT very rarely runs the whole span from chip 0 to the truncation 
chip. In actuality, for these TSPRT values of « and £, we found that with no modulation 
the expected value of ASN for the scheme 1s 95 chips. So, based on this, we would in- 
tuitively guess that modulation will cause the most damage if it occurs at the midpoint 
between chip 0 and chip 95. This 15, in fact, what we observe. 

Figure 46 on page 79 shows the ASN vs. modulation boundary under H, and 
Figure 47 on page 79 shows the PFA vs. modulation boundary. Note that, as with the 
FSS test, modulation has little effect on receiver performance under the H, hypothesis. 

We now summarize our TSPRT results for « =0.04, 6 =0.14. Figure 48 on 
page 80 shows the test power vs. |j - y | with and without modulation. Without mod- 
ulation, E[PD] is 0.92. With modulation, E[PD] drops to 0.89, but this is still greater 
than our desired value for the detection probability (0.86). This surprising result seems 
to indicate that out TSPRT 1s both overdesigned and robust. We only “request” that the 
test perform with E[PD] =0.86. It actually performs with E[PD] =0.92. If we then 
introduce modulation, the actual detection probability is decreased, but not enough such 
that its value is driven below our desired detection probability. This result can not be 
overemphasized: our original TSPRT scheme performs in accordance with our design 
values of o and fi even in the presence of data modulation. Figure 49 on page 80 sum- 
marizes the expected value of ASN vs. |j 4- y | with and without modulation. The effect 
of modulation on ASN is extremely small. 

The above results, while certainly encouraging, apply only to the particular case 
of œ =0.04, 8 =0.14. We wish to know just how robust and overdesigned the scheme 
1s for other values of « and f. 

We will pick a very stringent value of « and examine the TSPRT performance 
for various values off. Specifically, we let « — 0.01 and let f take on values from 0.01 
to 0.12. (Equivalently, our desired probability of detection ranges from 0.88 to 0.99.) 
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Figure 46. TSPRT ASN vs Modulation Chip Boundary: H, 
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Figure 47. TSPRT PFA vs Modulation Chip Boundary 
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Figure 48. TSPRT Power vs |; - y | With Modulation 
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Figure 49, TSPRT ASN vs |j+y| With Modulation 
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Figure 50 on page 82 shows the actual detection probability (from simulation) vs. the 
desired detection probability with and without modulation. Also shown is the optimal 
curve for the actual vs. design probability of detection. Note that even with this very 
Stringent value of =, the actual detection probability with modulation will satisfy our 
specification if our desired detection probability is less than 0.87. From Figure 50 we 
see that the degree to which the TSPRT 1s overdesigned increases as we relax the design 
value of detection probability. Figure 51 on page 82 displays the actual ASN vs. the 
design detection probability with and without modulation. We see that modulation has 
a negligible effect on the TSPRT ASN. 

In summary, modulation does not greatly affect out TSPRT sequential detection 
scheme. With a very strict desired false alarm probabilitv of 0.01 and a desired detection 
probability of 0.99, the actual detection probability in the presence of data modulation 
only drops to 0.95. If the desired detection probability is less than 0.88, the receiver will 
still perform in accordance with specifications. Note that this was not the case for the 
FSS test. For the FSS test, our results show that given any « and fj, modulation will 
drive the receiver out of specification. (The degree to which the FSS receiver is out of 


specification does, however, depend on the values of « and f.) 


D. THRESHOLD ADJUSTMENT 

We have seen that modulation causes no deterioration of the receiver performance 
under Z7, but does cause some degradation under Z,, particularlv for the FSS test. Since 
we know the exact form of the test statistic's pdf (given by (2.8) for H, and (4.7) or (4.9) 
for H,) we can redesign our decision processor to account for the modulation. We 
Epnsider the FSS test. 


Using the Q function, we can write the FSS false alarm and miss probabilities as 


Ра = а = (= | ze ~ 
См OM 
Js [v 
P miss = пр (928. =) 


n—1 
¿=0 
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Figure 50.  TSPRT PD With Modulation 
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Figure 51. TSPRT ASN With Modulation 
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We attempt to iteratively solve these equations to determine the values of t’ and the 
termination point M. Unfortunately, we find that for many values of « and f the 
equations have no joint solution. 

To investigate this dilemma, we first examine the test statistic pdf under H, (with 
modulation), given by (4.9), for various values of n, the number of chips in the inte- 
gration process. We have two unknowns to determine: :'/o?, (or z' in our normalized 
notation) and M. Equation (4.13) determines the value of z', since z' is the only un- 
known in (4.13). This value of z' is used in (4.14) to determine M. We try consecutive 
values of M in (4.14) until we find the first value of M that satisfies the equation. 

We find that sometimes no value of M exists which satisfies (4.14) once z' is deter- 
mined from (4.13). Figure 52 on page 84 shows a plot of (4.9), the test statistic pdf 
under H,, for various values of n. Note that the value of the density function at z 2 0 
increases as и increases. 

A particular example will illustrate the problem. Suppose we choose a — р = 0.01 
and we desire the solution of the pair of equations (4.13) and (4.14). Solving (4.13) we 
find that z’= 10.0777. Using this value of z' in (4.14), we attempt to adjust n until the 
equation 15 satisfied, but we find no solution. Instead of solving (4.14), let us try to de- 
termine a graphical solution. 

Noting that (4.14) is the cumulative distribution function for the pdf under H,, we 
can rephrase (4.14) in words: determine a value of n such that the area under the prob- 
ability density. function, given. by (4.9) and shown in Figure 52, from z=0 to 
2 = 10.077 15 less than or equal to 0.01. It turns out that the area from z=0 to 
z = 10.077 can not be made less than 0.01. We determined the exact value of the area 
Merete pdf between z — 0 апа 2 = 10.077 [ог а! values of n from n=0 to m= 1000. 
The results are shown in Figure 53 on page 85. We find that the minimum area ever 
possible between z 2 0 and z — 10.077 1s 0.0777 square units, which occurs when 
n — 280. Returning to (4.14), we conclude that if we attempt to adjust the thresholds, 
the best miss probability we can theoretically achieve (when х = 0.01) 15 5 = 0.0777, апа 
this occurs if we run the FSS for 280 chips. 


We draw two conclusions from our specific example. 


]. If we choose « =0.01 then, no matter how we attempt to adjust the FSS scheme, 
we can not theoretically achieve a value of p less than 0.0777. If we choose a value 
ог В greater than 0.0777, say 0.09 or 0.15, we should be able to adjust the FSS 
scheme such that it performs within specification when modulation is present. 


2. The value of n obtained for the smallest fj, n= 280 in our case, 1s the “optimal” п 
we can choose if we desire а В less than 0.0777. In other words, if we desire 
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Figure 52. Density Function With Modulation: H, 


В = 0.05 then, if we pick a test length of 280 chips, although we will not achieve 
our desired fj, any other value of M (less than or greater than 280) will perform 
even worse. 

Note that the example we chose above is particularly strict. For many values of a 
and f (other than the ones we used above) the two equations (4.13) and (4.14) will have 
a joint solution which gives us a successful adjusted FSS test. Only for very small values 
of B do the equations not have a solution. 

Figure 54 on page 86 shows the value of M, the adjusted FSS test length, for 
0.01 < а, x 0.20 with modulation. The optimal test length is found by solving (4.13) 
and (4.14) simultaneously for each pair (a,f). The flat portion of the curve, for small 
values of design miss probability, represent "best" choices for M, i.e., choices for M 


which will not satisfy our desired a and f but are, nevertheless, the best we can do. 
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Figure 55 on page 87 displays the required test length vs. desired miss probability for 
various fixed probabilities of false alarm. We expect the test to meet the specifications 
(x and f) for all values of the design miss probability except for points on the flat 
portions of the curves in Figure 55. For the flat portions, using the indicated value of 
M for the FSS test length will give us the best miss probability that can possibly be 
achieved, although we expect f) will not be satisfied. 

Figure 56 on page 88 shows the results of actual miss probability vs « and f with 
modulation present and thresholds suitably adjusted. Compare these results to those 
obtained without threshold adjustment, shown in Figure 40. We see that our threshold 
adjustment does indeed work where we expect it to, and does indeed fail where we pre- 
dicted it would. Prior to threshold adjustment, we see from Figure 40 that the test al- 
ways failed. Figure 57 on page 89 shows the actual miss probability vs. a for several 
fixed values of f. Note that if the desired miss probability is 0.1 or 0.2, the adjusted FSS 
test satisfies the performance criteria. If the desired miss probability is 0.01, we see that 
the test fails, but the curve in Figure 57 shows the best that is achievable. 

Looking at the results another way, Figure 58 on page 89 displays the actual miss 


probability vs. the design miss probability for several fixed values of a. For each value 
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оЁ х, the adjusted FSS test performs within specification until the design miss probability 


goes below a certain value. 

Finally, Figure 59 on page 90 shows the actual false alarm probability vs. a and fj 
for the FSS test with the thresholds adjusted to account for the data modulated received 
signal. Note that the results are similar to those obtained without threshold adjustment 


(Figure 43). Again, modulation does not significantly affect the test performance under 


the H, hypothesis. 
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Figure $4. Adjusted FSS Test Length 
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Figure 56. Actual FSS Miss Probability: New Thresholds 
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Figure 57. 


Figure 58. 
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Figure 59. Actual FSS False Alarms Probability: New Thresholds 
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У. CONCLUSIONS 


In this thesis, we have studied the effects of fading and modulation on noncoherent 
m-sequence acquisition schemes. 

We have considered two schemes: the fixed sample size test (FSS) and the truncated 
sequential probability ratio test (TSPRT). Most prior studies have assumed that co- 
herent demodulation is available during the acquisition process. We have dropped this 
unrealistic assumption since, in practice, m-sequence acquisition is generally performed 
bv the receiver prior to carrier recovery. Additionally, most prior research has assumed 
that the samples used in the decision processor are independent and identically distrib- 
uted. We do not reset the integrators during the correlation process, resulting in an en- 
hanced effective signal to noise ratio. 

The effects of channel fading on the performance of the synchronization tests were 
investigated. We considered a Ricean fading channel characterized by r, where r is the 
ratio of the received signal power in the direct component to the received signal power 
in the diffuse component. The probability density function for the receiver's test statistic 
was derived for the Ricean fading channel and found to be noncentral Chi-squared. It 
was determined that fading theoreticallv reduces the received signal's SNR bv a factor 
of r/[1 4 r 4- 2,201]. 

To simulate fading we considered two simulation techniques. The first was a Monte 
Carlo algorithm derived from the strong law of large numbers. It was found that the 
Monte Carlo method did not provide a very tight confidence interval for a 50-run sim- 
ulation. We then applied a numerical integration method that provided answers which 
converged within 25 runs. 

We examined the effects of channel fading on the performance of the original deci- 
sion processor. For the FSS test we observed that fading had virtuallv no effect on the 
false alarm probability, but drastically reduced the detection probability. With very se- 
vere fading, a detection probability of 0.99 can be decreased by more than 25%. Asr 
increases, the detection probability approaches the original value without fading. For 
the TSPRT we also observed that the false alarm probability was essentially unaffected 
by the fading, while the detection probability was severely reduced. For the TSPRT, the 
ASN under H, also remains essentially the same, but the ASN under H, increases as the 


fading becomes more severe. 


9] 


We then considered performance of the tests designed with the knowledge of the 
fading. We adjusted the sample size and thresholds of the decision schemes in accord- 
ance with the new test statistic density function. For the FSS test the performance cri- 
teria were then satisfied but the required sample size increased to a very large number 
when the fading was severe. As an example, for a = 8 = 0.01 and an SNR of -10dB, the 
required FSS sample size increased from 258 (for no fading) to 361, 509, 6617 and 8944 
when fading is characterized bv r= 20, 10, 1 and 0, respectively. Since the H, condition 
is unaffected by fading and occurs for all but a few of the uncertainty phases, adjusting 
the FSS thresholds greatly increases the acquisition time. Our results suggest that it may 
not be worthwhile to adjust the decision processor to account for fading. 

When we adjusted the thresholds for the TSPRT test, a similar result was observed. 
The updated test performed in the fading channel in accordance with specifications, but 
at the expense of a much larger ASN under //,. As an example, for a = 6 = 0.01 and an 
SNR of -10dB, the ASN under //, increased from 158 to 215 and 308 for r= 20 and 10, 
respectively. To reduce receiver acquisition time it 1s critical that the ASN be minimized 
under H,. In adjusting the test thresholds, the price to pay for achieving the probability 
of detection 1s a large increase in ASN under //,. We conclude that it might be prefer- 
able to apply a test designed under no-fading conditions; we may simply accept any de- 
gradation in detection probability caused bv the presence of fading in the channel. 

We then investigated the effects of data modulation on the performance of the syn- 
chronization tests. Again, an expression for the density function of the receiver test 
Statistic was derived. To see the effects of data modulation on the performance of the 
original decision processor, we first considered tests designed under no modulation and 
evaluated their performance in the presence of data modulation. For each scheme the 
modulation effect was found to be not severe. For example, И а = f/ 2 0.01 and SNR 
= -10dB, modulation causes almost the same effect on the FSS test as fading with r = 10. 
Modulation has virtually no effect on the FSS test false alarm probability, but does 
moderately reduce the detection probability. The degradation in the detection proba- 
bility tends to lessen as we relax the design criteria (a and f). 

We found that modulation had even less of an effect on the TSPRT scheme. For a 
stringent desired false alarm probability of 0.01, the scheme will meet any desired de- 
tection probability less than 0.87 even if data modulation is present. If the desired de- 
tection probability is greater than 0.87, the scheme will perform slightly worse than 
desired. For example, if the desired probability of detection for the original receiver is 
0.99, it is only reduced to 0.95 if data modulation is present. The ASN of the TSPRT 


is essentially unchanged with data modulation. We conclude that the TSPRT is rather 
robust in the presence of data modulation and threshold adjustment may be not war- 
ranted. 

We finally considered performance of the FSS test designed with the knowledge of 
data modulation. We found that there are values of (a,f) for which it is impossible for 
the adjusted FSS receiver to perform within specification. For these points we deter- 
mined the best sample size which resulted in the largest possible detection probability. 


For the other values of (a,f), the test can be adjusted to perform satisfactorilv at the 


expense of increasing the sample size. 
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APPENDIX GENERATION OF RICEAN RANDOM VARIATES 


To generate random variates from an underlying probability distribution, one starts 
with uniformly distributed random variates (on the interval [0,1]) and uses a suitable 
transformation to arrive at the desired distribution. Uniform variates are always the 
logical starting point: the theory 1s very well developed and many algorithms have been 
Written and thoroughly tested. 

Obviously computer-generated "random" numbers are not random; thev are, in fact, 
totally predictable. Furthermore, their predictability is often a desirable property [Ref. 
25 : pp. 132-145]. The best we can achieve with a computer random number generator 
is a sequence of numbers which appear to be random when subjected to standard statis- 
ucal tests. Generallv, the best random number generators put out sequences of numbers 
which have long periods, 1.e., which take a long time to start repeating. 

Making use of abstract algebra, random number generators of the form x, = 
f(x) have been developed which have the maximum possible period prior to repeating. 
The most commonly used maximum length generator 1s the congruential generator which 


generates uniform random variates u, by using the algorithm [Ref. 26: pp. 20-26]: 


ха = (ах; + с) (тойт) (4.1) 


(4.2) 


То get the period equal to m, we need to follow three rules [Ref. 22 : p. 15] 
l. c must be relatively prime to m 
2. a= | (mod p) for every prime p which divides m 


3. a z m (mod 4) if m 1s a multiple of 4 


To optimize the statistical properties of the random variates, we adhere to a fourth 
rule [Ref. 27 : pp. 155-156] 


4. the multiplier a should be larger than /m, preferably larger than m/100, but 
smaller than m — Jm 


The simulation programs were designed to run on a VAX computer using a base of 
2 and a word size of 31. Thus theoretically, the best choice for m (to result in the longest 


sequence period) would be m — 2". However we can not choose m= 2* because the 
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largest number generated by the algorithm (A.1)(A.2) 1s a(m — Il) - c. (This occurs 
when ax, +c=m-—1. In such a case (ax, + clímod m) will equal m — 1 and then x, 
will cause a number a(m— 1) + сто be generated.) If we choose m = 2* then, during the 
algorithm, at some point, the computer will want to work with a number larger than 23! 
(potentially as large as a(2?! —1) 4 c ), which will cause overflow. 

Rule 4. above suggests a possible choice for the constants a, m and c. If we choosc 
а ~ Јт then the largest number generated by the algorithm will be a(m— 1) 4 c 


~ mm -— m, We set this equal to the largest number obtainable in the computer: 
poo: (4.3) 


So, to the nearest power of two, a good choice for m is thus m — 2? . The reason for 
choosing a power of two 1s to facilitate choosing the constant c which, by the third rule, 
must be relatively prime to m. Any odd number can now be chosen for c. 

Summarizing the above results, we can satisfy rules 1.-4. and have a maximum ре- 
riod for our VAX computer by choosing: 


l. m = 2% = 1048576. This will be the number of uniform random variates gener- 
ated prior to repeating. 


2. Since we need a > „/т we must choose a > 1024. To satisfy a = | (mod2) and 
а = 1 (тоа 4) же сһооѕе а = 1029. 


3. For c pick any odd number, subject to 4. below. 


4. The maximum value generated bv the algorithm is a(2? — 1) + c= 2” + c. Since 
this number must be kept less than 2?!, we require that the number c be kept less 
than about one billion. 

For the remainder of this discussion, suppose we have written two separate uniform 
random number generators (e.g., by applying the choices outlined above but with two 
different values of c) which output independent variates U, and U,. We now use these 
to generate other needed distributions. 

It can be shown [Ref. 26 : pp. 67-68] that using the random variable U,, we can 
generate a random variate with an exponential distribution of parameter one by the 


transformation 


| 
E, = i к. | (4.4) 
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Now, with this exponential variate E, and the independent uniform random variable 
U,, we can generate two independent normal variates with distribution N(0,1) by using 


the procedure [Ref. 28] 


№ = J2E, cos2nzU, (А.5а) 
ХМ; = УБ нат (4.56) 


Suppose that instead of standard normal vanates N, and N, we desire normal 
variates described as N(s,o?) and N(0,o?). That is, we desire that both of our normal 
variates be independent but with variance o? (instead of one) and with means s and 0 
respectively. Such normal variates can be generated by using the basic properties of 


mean and variance [Ref 15 : pp. 96-103] as 
№ = oN + 5 (4.6a) 
Мв = 0“, (4.65) 


We have now generated two normal variates with №, ~ N(s,o?) and Ng ~ N(0,o?). 
To generate a Ricean distributed random variable, we need only to calculate [Ref | : pp. 


30-31] 


у = уха + ХК: (АЛ) 


Such a random variable will have a Ricean density function given by (3.2) and repeated 


below: 


М pg — 2199? Vs 
G о 
Recall from Chapter 111 that fading is usually characterized by specifying the value 
of r, where r is the ratio of the power in the direct component to the power in the diffuse 


component: 





with the constraint that total power has a normalized value of 1 (i.e., s? -- 2e? — 1). 
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We consider an example. Suppose that we wanted to generate independent variates 
from a Ricean distribution characterized by r= 1. One way to choose s? and o? would 


be 
52 =102 , 0 = 14 


so that s + 20? + |. So starting with our standard normal variates N, and N,, we 
calculate, using (A.6): 


: Ny 1 
Na а 


| Np 
D 


and then determine y from (A.7). 


927] 
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